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SUMMARY 
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The multiple quantum well (MQW) structure was first proposed in 1980 as a 
method to enhance the electron-hole ionization ratio in photodiodes beyond that typically 
found in bulk materials. Later in 1982, the doped MQW was introduced in order to further 
improve the electron ionization rate over that of holes. The desire for high performance 
optical detectors has resulted in several proposed MQW structures using different material 
systems in an attempt to optimize their gain, noise, and bandwidth characteristics. 

In this work, a detailed experimental investigation and analysis were performed of 
the physical properties of advanced semiconductor junctions. The analysis includes a study 
of (1) the difference in the structure-induced multiplication gain between doped 
GaAs/AJGaAs MQW and PIN junctions, (2) the effect of variations in the doping profiles 
on the properties of doped MQW structures, (3) the effect of surface treatments on the 
dark current and gain characteristics, (4) the spectral and time response limitations of the 
structures, and (5) a complete modeling of the junction physics for the different structure 
types All of these investigations were performed by conducting experimental 
measurements and theoretical simulations on new avalanche photodiode (APD) structures 
with built-in intrinsic (PIN), doped and undoped MQW structures. 

In order to properly interpret the experimental data provided by the measurements, 
one needs to understand the correlation between such data and the physical parameters 
used in designing the structures. To accurately determine such a correlation usually 
requires the ability to grow and fabricate a large sample of structures produced under 
very similar conditions. Unfortunately, this is not very practical or even possible during 
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material growth and fabrication. Therefore, developing theoretical models which 
accurately predict the relationships between the input and the output parameters is 
essentia] to understanding the physics behind the data. 

Most of the experimental data will be presented and analyzed for the first time in 
the GaAs/AlGaAs material system. The experimental results were compared to theoretical 
models, and were used to demonstrate, for the first time, the impact of the doping 
imbalance throughout the structure on the optical and electrical characteristics of a doped 
MQW structure. These models accurately predicted most of the external behavior 
displayed by these structures during experimental testing. In addition, various surface 
treatment techniques which enabled a dramatic reduction in the reverse bias dark current 
by as much as a factor of 1000 will be discussed. Furthermore, a new technique will be 
presented for improving the quantum efficiencies of these structures, and its effectiveness 
was verified through theoretical models. 


CHAPTER I 


INTRODUCTION 
1.1 Introduction 

Extensive studies have been performed on MQW junctions and structures because 
of their potential applications as avalanche photodetectors in optical communications and 

imaging systems. ' The role of the avalanche photodiode is to provide for the conversion of 
an optical signal into charge. Knowledge of junction physics, and the various carrier 
generation/recombination mechanisms, is crucial for effectively optimizing the conversion 
process and increasing the structure’s quantum efficiency. In addition, the recent interest 
in the use of APDs in imaging systems has necessitated the development of semiconductor 
junctions with low dark currents and high gains for low light applications. Because of the 
high frame rate and high pixel density requirements in new imaging applications, it is 
necessary to provide some front-end gain in the imager to allow operation under 
reasonable light conditions. Understanding the electron/hole impact ionization process, as 
well as diffusion and surface leakage effects, is needed to help maintain low dark currents 
and high gains for such applications. In addition, the APD must be capable of operating 
with low power, and low noise. Knowledge of the effects of various doping configurations 
and electric field profiles, as well as the excess noise resulting from the avalanche process, 
are needed to help maintain low operating bias and minimize the noise output. 
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To understand and quantify the physics of the different junctions, it is necessary to 
develop the structures into devices that could be studied on the macroscopic level and 
under different testing conditions. The various APD structures that have been developed 
for such a purpose include volume-doped wells/barriers multiple quantum well (MQW), 
delta-doped MQW, conventional PIN structures, and doped i-region PIN APD's. One way 
to quantify the gain and noise properties of an APD involves the measurement of the ratio 
of the impact ionization coefficients of electrons and holes. The ionization coefficient 
(measured in units of cm' 1 ) represents the inverse of the average distance traveled by a 
carrier between two consecutive ionizing collisions. In the majority of compound 
semiconductor materials, the electron ionization coefficient (a) is comparable to that of 
holes (P) and are both greater than zero. An ideal device where the hole-ionization 
coefficient (P) is equal to zero would have no multiplication noise and a performance 
similar to that of a photomultiplier tube (PMT) (gain ~ 10 s ) where multiplication takes 
place with just one type of carriers. In bulk GaAs materials, k (=a/p) is equal to about 
1 .67. In bulk Silicon, k approaches a value of 20, but still does not compare to the single- 
carrier multiplication characteristics of PMT. Despite the advantages of PMTs, these 
photodetectors tend to be cumbersome, have low quantum efficiencies, and operate at 
very high voltages (1,000-10,000 V). On the other hand, APDs are small, rugged, have 
relatively low-cost, and low operating voltages (10-100 V), as well as high quantum 
efficiencies. The desire for high performance optical detectors has resulted in several 
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proposed MQW APD structures using different material systems in an attempt to optimize 
their gain, noise, and bandwidth characteristics. 

In theory, the value of k in a MQW structure is much higher than that typically 
found in bulk materials since it is possible to design a structure which selectively enhances 
the energy of one carrier over another. In order to design such a structure, as well as 
optimize the other properties of a photodetector such as gain, quantum efficiency and time 
response, one needs to understand the various physical mechanisms that take place inside 
the semiconductor structures In addition, the impact of various growth conditions on the 
structure’s behavior must be understood in order to consistently grow devices with similar 
optical and electrical properties. 

Different experimental techniques including current-voltage (IV), capacitance- 
voltage (CV), gain, excess noise, spectral and pulse response measurements were used to 
extract information on the various optical and energy band related processes and 
mechanisms that control the performance of these structures. Gain studies, for example, 
provide insights into carrier multiplication via impact excitation of confined electrons out 
of the narrow-gap semiconductor layers. Similarly, noise studies allow the measure of the 
ratio of the impact ionization ratio of electrons to that of holes, which is determined by 
various structure parameters such as welL'barrier widths and the doping concentrations in 
the wells. CV data is used to determine the carrier concentration profile versus depth 
throughout the structure, while spectral response provides insights about the structures’ 
quantum efficiency and its sensitivity to various electromagnetic wavelengths. Finally, time 
response and device bandwidth measurements provide information on the carrier transit 
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time in the electric field region, carrier diffusion time in the undepleted absorption region, 
as well as carrier trapping at the heterojunction interfaces. 

1 .2 APDs in Receivers 

One common application of the avalanche photodiode structure involves using the 
device as an optical receiver in communication systems. The increasing requirements of 
high bit rate fiber communication networks require the development of suitable high 
performance optoelectronic components. The photodetector (APD) is used at the 
receiving end of the fiber to convert the optical signal into an electrical current which can 
be processed to recover the original data. For low bit rate and short range applications 
such as a local area network (LAN), the device performance requirements are only 
moderately stringent. For high bit rate, long distance applications and bulk data 
transmission (high resolution graphics/sound, high definition television), the device 
performance requirements become far more demanding in terms of speed, size, reliability, 
etc. 

High speed electronics have made it possible to encode/decode data at 
continuously faster rates. The recent development of state of the art narrow linewidth 
lasers and low dispersion optical fibers has made it possible to transmit light at high bit 
rates for longer distances with minimum degradation in signal quality. With these 
developments, the limiting factors in fiber optical systems seems to be determined by the 
receiver sensitivity and speed performance. This provides considerable incentive to study 
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and develop high performance optical detectors that are designed to handle such stringent 
requirements. 

1.3 APDs in Imaging Systems 

Another related application of APDs is in the development of imaging systems to 
be used both in the visible and non-visible portions of the spectrum. While fiber optic 
communications are usually limited to the 1.3-1.55 pm spectral range, imaging systems 
can be developed using a wide range of materials to operate in a variety of optical ranges 
depending on the application in question. The role of the APD in such a system is to 
capture the optical images and to amplify the signal using its built-in gain properties. One 
such system is currently under development by the Georgia Tech Research Institute to be 
used in a high definition television (HDTV) imaging system. This system uses a 1920 x 
1080 imaging array of superlattice GaAs/AlGaAs APDs as an image capture mechanism 
operating at wavelengths below 1 pm. 

1.4 Background and History of Semiconductor Junctions used in APDs 
1.4,1 Common Semiconductor Materials 

Silicon APDs are one of the most widely used photodetector devices in the 
spectral range below 1.1 pm. Extensive research has been conducted in this area and the 
behavior of Si junctions is fairly well understood 2,3,4 . Unfortunately, these devices are not 
very useful for fiber optic communication because of the high signal dispersion and 
attenuation at these wavelengths. Dispersion in fibers is zero around 1.3 pm, and 
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attenuation is minimized near 1.5 pm thereby requiring the use of semiconductor materials 
that will operate more efficiently at these wavelengths. Ge APDs have been developed to 
operate at these wavelengths, but they were found to have high dark currents and equal 
electron and hole impact ionization rates which limits their gain/noise properties. For that 
reason, recent work has been focused largely on lattice-matched III-V materials with 
adjustable energy gaps that can be customized to optimally operate at a wide range of 
wavelengths throughout the visible and near infrared part of the spectrum. Such material 
systems include ternary alloys such as Al x Gai. x As and quaternary alloys such as In x Gai. 
x AsyPi- y where x and y are the mole fractions of the group III and V elements 5 . The 
complex MQW structures require a great degree of doping control as well as high 
interface quality during material growth. The progress realized in GaAs/AlGaAs growth 
by molecular beam epitaxy (MBE) during the past fifteen years makes this material system 
an ideal candidate for MQW structures. 

1.4,2 Comparison of Semiconductor Junctions used in APDs 
1 .4.2. 1 PN/PIN Junctions 

Some of the first semiconductor junctions that were developed to be used as 
photodetectors were the PN/PIN structures. The PIN “APD” was first proposed by Read 6 
in 1958 and it operates in a manner similar to a PIN “photodiode”, with the exception that 
the primary photocurrent caused by the photoelectric effect is amplified within the same 
device. An APD-based receiver is more sensitive than a PIN photodiode-based receiver, 
provided the APD has sufficient bandwidth for the application 7 . However, an APD’s 


6 



bandwidth is significantly less than that of an equivalent PIN photodiode In addition, as 
previously mentioned, there is excess noise associated with the avalanche multiplication 
process. This can limit the maximum useful gain of the APD to the point where the APD 
noise becomes comparable to the first stage amplifier noise 7 . In addition, the presence of 
high dark currents usually presents a limiting factor preventing the further increase in a 
device’s photocurrent gain beyond avalanche breakdown. The PN/PIN junctions are 
simple to analyze using a one-dimensional approximation which simplifies the equations of 
state to a single spatial variable and enables closed-form solutions of the differential 
equations. A typical PIN structure is shown in Figure 1-1, and is composed of a highly 
doped p + layer, followed by an intrinsic layer and a highly doped n + layer. Electron or hole 
injection is accomplished by shining light on the p + or n + regions, respectively. Single- 
carrier injection can be obtained by designing the thickness of the contact region to be 
small enough to allow diffusion of carriers to the high field region, and yet large enough to 
prevent direct absorption of light in that region. In GaAs, the ideal contact region is 
generally between 1 and 3 pm depending upon the wavelength of the incident light. 8 
1.4 2.2 Doped/Undoped MOW Junctions 

The MQW structure (Figure 1-2) was first proposed in 1980 by Chin et al 9 and 
Capasso et al. 10 (1982) as a method to enhance the electron-hole ionization ratio beyond 
that in bulk materials. Later in 1982, Blauvelt" proposed the doped MQW structure which 
incorporates built-in field layers in the avalanche region in an attempt to further increase k. 
Since then, various new ideas have been proposed to help increase (or decrease) k. These 
proposals are all based on one or more of the following principles: 
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Figure 1-1: (a) Schematic of a PIN photodiode under reverse bias, (b) electric field profile 
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• The elimination of the “feed-back” process by confining the carriers in potential wells 
formed between two heterojunctions 12 . 

• The utilization of the difference in the ionization energies and the quasi-electric fields 
between the electrons and holes. These differences are incorporated into the material 
by gradually changing the energy gap and creating a “graded-gap APD” 13 . 

• The utilization of the asymmetries between the conduction and valence band offsets 
present in composite III-V semiconductors (i.e, MQW, “staircase” APDs, etc ) 14,15 ’ 16 

• The utilization of periodic doping layers in the MQW structure (PN heterojunction, 
doped-barrier MQW, doped-well MQW, delta-doping) 17 ' 18 

In 1982, Capasso conducted experimental measurements on a GaAs/AlcusGao ssAs 
MQW structure where he observed a value of a three times larger than that in bulk GaAs. 
He explained the results by suggesting that the conduction band discontinuity helps to 
contribute to the total ionization energy of electrons by reducing the impact ionization 
threshold energy by a value equal to AE C . On the other hand, the much smaller valence 
band discontinuity (AEv) does not significantly reduce the impact ionization threshold 
energy of holes. This results, according to Capasso, in increasing the value of a while 
keeping that of P relatively unchanged with respect to that found in bulk GaAs. However, 
as Aristin 8 pointed out, according to the conservation of energy principle, there is no 
reason for the electron ionization coefficient to be increased over any given period since 
the total gain of energy is zero. This has to be true since the energy gain at the 
AlGaAs/GaAs interface is offset by an equal loss at the GaAs/AJGaAs interface. 
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Therefore, from a conservation of energy perspective, it would seem as if an overall gain 
of energy is only possible in structures where the energy gap is continuously graded 13 . 

In 1987, Brennan et al. 19 were the first to explain the origin of the observable 
increase in the electron impact ionization coefficient a M Qw in the MQW structure. Using a 
simplified analytical expression of a(E), they have shown that, by replacing the square 
potential wells by a series of Fourier functions V(z), the final expression for <x M qw is 
always higher than that found in bulk GaAs materials. Such an increase is the result of the 
exponential dependence of a on the periodic electric field in the MQW structure. 
Therefore, even though the science itself is not yet very well understood, it became clear 
that the artificially enhanced ionization process inherent in MQW structures holds the key 
for designing optoelectronic devices with properties that could be optimized for a given 
application. 
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CHAPTER II 


THEORETICAL MODELING AND PHYSICS OF SEMICONDUCTOR JUNCTIONS 

2.1 Introduction 

The modeling of the various semiconductor junctions and APD devices was 
conducted using Silvaco’s AtlasII® device modeling framework which works by solving 
the basic Poisson and continuity equations for electrons and holes. Poisson’s equations 
describe variations in the electrostatic potential and how they relate to local charge 
densities while the continuity equations represent the effect of carrier transport and the 
recombination/generation processes on the material charge densities. 

2.2 Basic Semiconductor Equations 
2.2.1 Poisson and Continuity Equations 

The electrical properties of semiconductor junctions can be described by the 
following Poisson’s equation : 

div(eV y/) = -q{p - n + N» - N' A ) - p F [2-1 ] 

as well as the following carrier continuity equations for both electrons and holes : 
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= -VJ r -U, 
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where y/ is the intrinsic Fermi potential, N* D and N A represent the donor and acceptor 
ionized impurity concentrations, p F is the fixed charge density, U n and U p are the net 

electron and hole recombination rates, respectively. 

The modeling program solves the above three partial differential equations for the 
electrostatic potential, y / , and for the electrons and hole concentrations, n and p, 
respectively. 


2.2.2 Position-dependent Current Density Equations 

The electron and hole current density equations can be written in terms of carrier 
concentrations, n and p, carrier mobilities, /j„and p p , and the quasi-Fermi potentials for 

electrons and holes, ^„and <p p : 


J n = 

J P = -MpPV</> P 


[2-3] 


where, 


*■ = 

i 

q E fp 


[2—4] 
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and the Fermi energies are expressed as: 


E m =E e + kT\njr + kT\ny n 
E r r=E,-kT\n-£- + kT\r,y r 


[2-5] 


The last terms in the above two equations are due to the influence of Fermi-Dirac statistics 
and are given later by equations [2-18] and [2-19] in section 2.3.2 . In the case of 
Boltzman’s statistics (see section 2.3.1), y n = y p = 1. In multi-layered heterojunction 
structures such as the ones used in MQW APDs, the energy band diagram is not uniform. 
The position-dependent conduction and valence band energies can therefore be written as: 

E c = q{Vo-¥)~X 
E v = q(W 0 -V)~ X~E g 


where, 

X is the position-dependent electron affinity, E g is the position-dependent bandgap, and vj/ 0 
is some reference potential which can be selected in the form: 


Xr kT N„ Xr+E, kT N w 

A r , I CT __ ^ r & J ^ VI 


^o = — + — In 
q q n, 


[2-7] 


<7 q ”ir 

where n ir is the intrinsic carrier concentration of the arbitrarily selected reference material, 
and r is an index indicating that all of the parameters are taken from the reference material. 

Consequently, by combining equations [2-3]-[2-7] with [2-18] and [2-19] the 

following position-dependent drift-diffusion current densities can be obtained: 


14 


[ 2 - 8 ] 


J n - kT/u„'Vn-q/i n nV 


kT * kT N 
w + — In v + — + — In — 
V q ,n q q n„ j 

f 


kT y + E kT N 

J p =-kT Mp Vn-q MpP V I v lny,+ In — 

\ q q q n «J 


2.3 Carrier Statistics 

The electron and hole concentration densities can be defined using Fermi-Dirac 
distributions and a parabolic density of states giving : 


n=N c F V2 {±-[E Fn -E c }} 

k ! [2-9] 

P=N v F m {-[E y -E Fp ]} 


where N v and N c are the effective densities of states in the valence and conduction bands, 
E v and E c are the valence and conduction bands energies, and Ef„ = -qf n and E Fp =-pfp are 
the electron and hole Fermi energies. The Fermi-Dirac integral of order one-half is defined 
as . 






.1/2 


+ exp 




[2-10] 


2.3.1 Boltzman Statistics 

For the range of operation of most semiconductor devices, the electron and hole 
concentration equations can be simplified using Boltzmann statistics as follows . 
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[ 2 - 11 ] 


1 <7 

n*N c exp{— [E Fn -E c ]} = n, expM-(y/ - <j> n )) 
kT kT 

p*N v exp {jf[Ey ~ E Fp \) = n, exp[-^(^ - y/)] 

In the case when band-gap narrowing can be neglected, the intrinsic carrier concentration 
is expressed by : 


n,(T) = JN C N V exp (~E g / 2kT) [2-12] 

The band-gap and effective density of states have the following temperature dependencies 
according to Sze 20 : 


E g (T) = E g (0)- 


aT 2 
T + p 


(300) + a 


300 2 
300 + P 


T 2 

T + p 


N C (T) = 2 


2nm dt kT\ 

h 2 J 


3/2 



3/2 

N c (300) 


N V {T) 




n 3/2 


300/ 


AV(300) 


[2-13] 


[2-14] 


[2-15] 


mde and rridh are the density-of-state effective masses of the valence and conduction bands, 
respectively. 

The intrinsic Fermi potential is given by : 


~Wi 


Ec 



kT, N c 
— In — - 


2 Ny 


[2-16] 
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Under Boltzman statistics, the correlation between the mobilities and diffusivities in 


Equation [2-4] is given by the following Einstein’s relationships : 



[2-17] 


2.3.2 Fermi-Dirac Statistics 

The AtlasII modeling program uses both Boltzman and Fermi-Dirac statistics. The 
form of the density of states equations can be adjusted by introducing degeneracy factors 
y „ and y p given as : 


Yn - ^1/2 {^^[^Fn £ c ]}exp] 




y p- h \rA~j.[E r - £,,])exp 


]_ 

kT 


--L-{E v -E F p ) 


and the density of states equations can be rewritten as . 


n=N c y n ex p 


n = N c y „ exp 



[2-18] 


[2-19] 


[2-20] 


[ 2 - 21 ] 


When Fermi-Dirac statistics are used, Einstein’s relationships must be modified as follows 
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[ 2 - 22 ] 


2.3.3 Carrier Recombination Statistics 

The modeling program supports the following three recombination mechanisms : 
a) Shockley-Read-Hall: 


U srn ~ ' 


pn-n k 


In + n u exp| 


E,-E, 

kT 


+ T. 


p + n H ex 




-E, 

kT 


[2-23] 


b) Auger: U Augtr = c n ( pn 2 - im], ) + c p (np 2 - pn 2 t ) [2-24] 

c) Optical Recombination: 

U optica ,=C opl (pn-n l) [2-25] 

where Ej is the intrinsic Fermi energy, E, is the trap energy level, and n; e is the effective 
intrinsic concentration. Copt, Cn, and Cp are material recombination parameters for the 
GaAs/AlGaAs systems. The electron and hole lifetime parameters, r„ and r used in 
Equation [2-19] can be functions of the impurity concentrations as follows 21 : 
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[2-26] 


*A x ,y) 


no 


1 + N(x,y)/N SRH -„ 


i p (x,y) = 


1 + N(x,y)f Nsw-,, 


[2-27] 


where N(x,y) is the localized total impurity concentration. N S rh-«, Nsrh-p, t p0 , and r n0 are 
material parameters. 


2.3.4 Carrier Generation 

In addition to photogeneration of electron-hole pairs, the model can be modified to 
incorporate other carrier generation mechanisms such as impact ionization and generation 
due to band-to-band tunneling using the models presented in the following two sections. 

2.3.4. 1 Impact Ionization 

Avalanche multiplication due to impact ionization is the most important process 
during junction breakdown. Under the effect of an electric field, a carrier (electron or 
hole) acquires sufficient energy that, upon impact with the lattice, an electron-hole pair (e- 
h) is produced. These new carriers are accelerated by the field and can themselves acquire 
high-energy, causing more e-h pairs to be created. If insufficient energy is transferred to 
the lattice, impact ionization does not occur, and the energy is lost, usually by heating. In 
addition, carriers can lose energy through non-ionizing scattering events before acquiring 
sufficient energy to ionize. The scattering rate can be different for electrons and holes. 
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The electron-hole generation rate due to impact ionization in the field region was 


modeled according to Selberherr 22 as follows : 


G = c„^ + a 

q 1 


\J. 


[2-28] 


where a„and a p are the electron and hole ionization rates which can be expressed as a 
function of the electric field as : 


a 


n 



= a„ exp 



[2-29] 


a 


p 


= a, exp 



r E c nt ^ 

P,~ 


p 



l E j 



[2-30] 


where E is the component of the electric field in the current flow direction. All other 
impact ionization parameters are material dependent and are given in Appendix A for the 
AlGaAs material system, 


2. 3. 4. 2 Band-to-Band Tunneling 

In the presence of a high electric field, electrons can tunnel from the valence band 
to the conduction band in a reverse-biased junction. The criteria for tunneling are met in a 
high-field junction where the depletion width is small, and the potential barrier (in a MQW 
structure) is very thin. This process is shown for a p-n junction in Figure 2-1 where the 
conduction band edge on the n-side (E cn ) drops below the valence band edge on the p-side 
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(Evp), providing empty energy states for the electrons to tunnel into. Tunneling will 
increase the electron generation rate yielding larger reverse currents. The quantum 
mechanical tunneling transmission probability is given by 23 : 


T,= 


1 + 


El sinh 2 kW 
4 E(E 0 -E) 


-i-i 


[ 2 - 31 ] 


where Eo and W are the barrier height and thickness, E is the energy of the carrier, and: 


K & 


2 m(E 0 - E ) 


[ 2 - 32 ] 


The tunneling current density is given by 


24 . 


J, = 


Jim* q 3 £V 
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3 qfr 


[ 2 - 33 ] 


J 


where V is the applied voltage, m* is the effective mass, E g is the bandgap, and ^ is the 
junction electric field. The resulting tunneling generation rate can be expressed by the 
Klaasen model as. 


G hb , — ^bbt^ r ex P[ 


B, 


bbt 


[ 2 - 34 ] 


where E is the electric field magnitude, At*t, Bbbt, and y are material constants. 
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2.3.5 Carrier Mobility Modeling 


Throughout the junction, carriers are accelerated by the local electric field, but 
loose some momentum due to various scattering mechanisms. These scattering processes 
are a result of impurities, lattice vibrations (phonons), other carriers, surfaces, and other 
material imperfections. The macroscopic mobility parameters used in the current densities 
equation [2-3] express the effects of these microscopic phenomena. These mobilities are 
functions of the local electric field, the lattice temperature, the doping concentration, etc. 
At low enough fields, the mobility has a characteristic low-field value usually denoted by 
|ic. This is a result of the fact that carriers are almost in equilibrium with the lattice. The 
value of |io is impacted by phonon and by impurity scattering. When the temperature of the 
lattice increases, phonon scattering increases, and the value of po decreases. Similarly, 
when the doping concentration increases, impurity scattering increases which causes po to 
decrease as well. On the other hand for high fields, the carrier mobility declines. This is 
because carriers with high energies can take part in a wider range of scattering processes. 
The mobility models used in our simulation are both doping- and field-dependent. The 
analytical mobility functions and data are provided in section A. 5 of Appendix A. 

2.4 PN/PIN junction equilibrium electrostatics 
Figure 2-2 shows a one-dimensional step pn junction which will be used for the 

purpose of the analysis. In order to obtain a closed-form solution for such a junction, the 

following assumptions must be made 25 : 

a) The device is one-dimensional (see Figure 2-2) 
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b) At x=0, there is an abrupt metallurgical junction. 

c) The p and n layers are uniformly doped with a step junction from N A to N D (Figure 2-3) 

d) The p and n contacts are perfect ohmic contacts and are far separated from the 
junction. 

Under equilibrium conditions, the “hypothetical” p and n initial carrier 
concentrations can be depicted as shown in Figure 2-4. Because of the concentration 
imbalance between the two sides of the junction, the carriers would begin to diffuse in 
order to make their distribution more homogeneous throughout the material. Therefore, 
the holes will tend to diffuse from the p- to the n-side and the electrons from the n- to the 
p-side leaving behind ionized acceptor (N” A ) and donor (NT d ) atoms. As a result, a net 
charge density (see Figure 2-5(b)) will be created by the reduction of the majority carrier 
concentrations. 

The net charge density leads to the creation of an electric field and a built-in 
potential difference across the junction which inhibit further diffusion of the majority 
carriers. From Gauss’s law, the electric field can be determined as follows: 

£(*) = — jp(x)dx, (V/ cm) [2-35] 

K,e 0 i 


where, 

K,= relative semiconductor dielectric constant 

eo = 8.854 x 10* 14 , (farad/cm) 

p(x) = q(p - n + N D - N a ) , (coulombs/cm 3 ) 
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Figure 2-2: One-dimensional step junction 



Figure 2-3: Step junction 25 



(a) (b) 


Figure 2-4: Initial hole (a) and electron (b) carrier concentrations across the junction 





i.e, p(x) is the imbalance between the charge carriers and the ions. Therefore, the electric 
field can be obtained through a graphical integration of Figure 2-5(b) which results in the 
profile shown in Figure 2-5(c). Subsequently, the potential gradient within the depletion 
region and across the device can be calculated from electromagnetic field theory as 
follows : 

V(x) = -jf(x)dx [2-36] 


where the potential reference was chosen such that V(-oo)=0. The potential diagram is 
shown in Figure 2-5(d). Therefore, there exists a built-in potential Vbi across the depletion 
region at equilibrium which results in energy band bending between the two sides of the 
junction and an energy difference equal to qV b i as seen in Figure 2-6. 

The net charge density can be derived from the energy band diagram using: 


and. 



which yields, 


P = 


f K ’ £ °) 


v q ) 

1 dx 2 ) 


[2-37] 


[2-38] 


[2-39] 
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p-Bulk region 


Depletion region «-Bulk region 


Figure 2-6: Energy band diagram for a pn junction at thermal equilibrium 
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2.4.1 Built-in Potential Calculation 


In order to derive an expression for the built-in potential, Vbi, we will need to 
analyze the junction at thermal equilibrium with no applied bias and no net current flowing 
across the junction: 

dn 

J N = 0 = ^N\dhft + J N\diffiaion ~ q/jn£, + qD N — [2-40] 


which yields. 


£ = - 






n)\dx 


'kT 
v q) 


lV dn 


nJ\dx 


The built-in voltage can therefore be calculated from equation [2-18] 
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kT ” ( r* ) dn 


r an 
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with, 


which yields. 


” p = »(-*) = ^ I 
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K, = In! 


N d N a 

n 2 


[2-41] 


[2-42] 

[2-43] 

[2-44] 


[2-45] 


2.4.2 Depletion Region Calculation 

In general, Poisson’s equation [2-1] is not easily solved in “closed form” for most 
junctions because p and n are in turn functions of the unknowns variables V and x. In 
order to obtain an explicit solution of V as a function of x, one must make a set of 
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assumptions which is referred to as the depletion approximation. This approximation 
assumes that the mobile carrier concentrations (n and p) are small compared to the donor 
and acceptor levels in the depletion region, and that charge neutrality exists elsewhere in 
the region : 

1 . N a » np or p p , i.e, p = -qN A for -Xp < x < 0. 

2. N d » n„ or p n , i.e, p = qN D for 0 < x < x„. 

3 . r = 0 for x > x„ and x < -x p 

The depletion approximation reduces Poisson’s equation to : 



d£ gN D 

— = * D for 0 < x < x 
dx K s e 0 

[2—46] 

and, 

dt, qN D 
dx = K s e q 

[2-47] 


The above equations can be solved for the electric field keeping in mind that the field is 
zero in the bulk regions and at the edges of the depletion region. This results in the 
following p- and n-side depletion region approximations of the electric field: 


£(*) = ( X P + *)» for -x p <x< 0 [2—48] 

and, 

= for 0<x<x n [2-49] 

Since the electric field must be continuous at x=0, the above two equations can be reduced 
to . 

H A * P = H D x n [2-50] 
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which states that the areas in Figure 2-5(b) are equal and that the total negative charge 
must equal the total positive charge throughout the junction. Using Equation [2-36], we 
can derive the expressions for the potentials at both sides of the junction which are given 
by. 

V(x) = ^~(x,+xY for - x <x< 0 [2-51] 

-aN 

and, V(x) = - jr^(x n - x) 1 +V bj for 0<x<x n [2-52] 

2 K,e 0 

The depletion layer width can now be calculated by making use of the boundary condition 
for the potential function, i.e, V(0')=V(0 + ), and with the aid of Equation [2-50], we have: 


and, 


2K s e 0 V b , N a 

q N d (N a+ N d ) 

2K s e 0 V bl N d 

q N A (N A+ N D ) 


W = x+x n = 

P n 


2 K,e 0 V bl (N a +N d ) 


il/2 


N a N d 


[2-53] 

[2-54] 

[2-55] 


which is valid for a two-sided abrupt junction. In the case of a one-sided junction (p + n or 
pn + ), Equation [2-55] reduces to: 


W = 


2 W*. 
B 


[2-56] 


where N b =Nd or N A depending on whether N A »N D or vice versa. The previous equations 
give the depletion width under built-in bias. In the case of externally applied voltage, V b j is 
replaced by (V b j-V A ) where V A is the applied bias in volts and is positive for forward bias 
and negative for reverse bias. 
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2.4.3 Junction Capacitance 


The junction capacitance of the structure can be approximated by that of a parallel 
plate capacitance given by: 

[2-57, 

where Cj is a function of the applied bias and is decreased under reverse bias due to the 
increase in the value of the depletion width W. The above relationship can be used to 
accurately determine the depletion width of a junction with a known capacitance. 
Experimentally, the capacitance of a junction can be measured as a function of bias using a 
capacitance meter such as the HP4277A LCZ unit that was used in our experiments. 
Equation [2-56] can then be used to determine the doping profile N B as a function of 
depletion width W assuming a one-sided depletion approximation. 


2.4.4 Multiplication and Ionization Coefficients 

As was mentioned in Chapter I, the electric field profile in a pin photodiode is 
constant throughout the intrinsic region. The values for the electron and hole impact 
ionization coeficients, a and P, can be calculated using the multiplication parameters as 
follows: 


a(E) 


i K(n-i 


In 


r M n (V)^ 


wM n (V)-M(V) \M p (V)J 


[2-58] 


/?(£) = 


1 


In 


r M p {V)^ 


wM(V)-M n (V) 


[2-59] 


where M„and M p are the electron and hole multiplication gains and are given by: 
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M(V) = 


[2-60] 


i»(v)-i D <n 

^ phO 

where, 

Iph(V) = photocurrent at the applied voltage V 
Id(V) = dark current at the applied voltage V 
Ipho = photocurrent at no applied bias (unity gain) 

2.5 Numerical Methods 

Based on the previously described theories, the modeling program uses different 
numerical methods for calculating the solutions to semiconductor device problems. Device 
operation is simulated using a set of anywhere from one to six coupled, non-linear, partial 
differential equations. The program produces numerical solutions to these equations by 
calculating the values of unknowns on a series of mesh points within the device structure. 
The original continuous model is converted to a discrete non-linear algebraic system that 
behaves approximately in the same manner. The non-linear algebraic system of equations 
is solved using an iterative procedure that refines consecutive estimates of the original 
guess. The iterative process continues until each successive correction is small enough to 
meet convergence criteria, or until it becomes apparent that the procedure is not going to 
converge. 

For each model, there are basically three types of solution techniques: (a) de- 
coupled (Gummel), (b) fully coupled (Newton) and (c) Block The decoupled technique 
such as the Gummel method will solve for each unknown in turn while keeping the other 
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variables constant. The fully coupled techniques such as the Newton’s method solve the 
total system of unknowns together. The Block method is a combination of the two. In 
other words, it will solve some equations fully coupled, while others are de-coupled. In 
general, the Gummel method is useful when the system of equations has linear 
convergence and is weakly coupled. The Newton method is useful in the case of quadratic 
convergence and when the system of equations is strongly coupled. The Block method has 
the advantage of faster simulation time over that of Newton, but is not as good as the 
Gummel in providing initial guesses to the solutions. 

A good initial guess for the variables to be evaluated is crucial for obtaining 
convergence. When no previous solutions exist, the initial guess is usually calculated from 
the supplied structure parameters. For example, the initial guess (at zero bias) for potential 
and carrier concentrations can be made using the specified doping profile, etc. During bias 
ramps such as the ones used in our IV and CV calculations, the initial guess for any bias 
point is obtained by an extrapolation of the two previous results. The extrapolation 
method for the initial guess will generally give good results when the variables measured 
(such as the current in IV) have linear characteristics. However, problems may arise when 
the variable’s behavior becomes highly nonlinear or change quickly such as is the case near 
breakdown or threshold. This will normally require repeated simulations to determine the 
threshold point at which the voltage steps must be reduced to allow proper convergence . 26 
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3.2 Simulation results for a doped/undoped MQW structure 



Figure 3-1: Electric field profile for an undoped 9-wells MQW APD under no bias 
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Figure 3-2: Conduction/valence band diagrams for an unbiased, undoped MQW 
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Figure 3-3: Difference in built-in e-field profile between a PIN and an undoped 

MQW 
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3.3 Simulation results for an undoped MOW structure 


Figures 3-6 and 3-7 shows the simulated electric field profile and the energy band 
diagram for a 9- well/ 10-barrier undoped MQW APD. The cap and bottom GaAs layers 
were 0.85 pm thick and were doped at lxlO 1 * cm* 3 . The Alo .42 Gao.jg As barriers and the 
GaAs wells were 500 A each and were separated from the top and bottom layers by 0.2 
pm thick intrinsic GaAs layers. The structure was designed to be symmetric with respect 
to a plane at 1.525 pm from either surface. This resulted in a symmetric electric field 
profile with a peak value at zero bias of around 45 kV/cm at the two doped/undoped 
GaAs interfaces. The field in the GaAs wells is uniform across the structure and has a 
value of about 9.3 kV/cm. The magnitude of the field in the AlGaAs barriers is lower at 
8.5 kV/cm as a result of the higher dielectric constant of the material as can be seen from 
Equation [A-7] in Appendix A. 

The diagram in Figure 3-7 shows the zero-bias valence/conduction band diagram 
of the structure. Since the built-in voltage is determined by N A , No, and n„ it is equal to 
about 1.375 V as is the case for the PIN structure. Since the built-in bias is equal to the 
spatial integral of the electric field across the junction, one would expect the electric field 
magnitude to be slightly higher in the GaAs wells than that in the intrinsic region of a PIN 
APD, and lower in the AlGaAs barriers. This is shown schematically in Figure 3-8. This 
behavior helps to explain why the avalanche breakdown voltage in a MQW APD tends to 
be slightly lower than that of a conventional PIN with the same dimensions. This is 
because the higher field magnitude in the GaAs wells tends to induce the impact 
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Figure 3-7: Conduction/valence band diagrams for an unbiased, undoped MQW APD 
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ionization process at a lower applied bias than would be possible if the MQW structure 
was not present. 

3.4 Simulation results for a doped MOW structure 

Figures 3-9 and 3-10 show the electric field and energy band diagrams for a 
doped-wells MQW structure at zero bias. The top p + and bottom n + layers were heavily 
doped at 3xl0 18 cm' 3 . All nine wells are doped with 50 A wide p and n layers at 1.5xl0 18 
cm' 3 . This creates localized high field regions throughout the structure with a peak value 
of around 95 kV/cm at zero bias as shown in Figure 3-9. 

Figure 3-10 shows the corresponding band diagram where the effect of well 
doping on the band structure is outlined. The additional energy drop per well as a result of 
this particular doping configuration is equal to about 0.06 eV. It is possible to further 
enhance the energy gain by increasing the doping density as well as the widths of the p and 
n doping layers. 

Figure 3-11 shows the field profile at a reverse bias of 4 V for the same doped-well 
APD previously described with the exception that the doping in the GaAs wells is such 
that p=1.8xl0 18 cm' 3 and n=1.5xl0 18 cm' 3 . The doping imbalance in the wells ruins the 
symmetry in the electric field profile and results in undepleted MQW stages throughout 
the structure as shown in the band diagram in Figure 3-12. The undepleted or “inactive” 
stages are low-field regions which present highly resistive barriers to the flow of electrons 
and holes across the device. Furthermore, the AlGaAs barriers tend to reduce the internal 
quantum efficiency of the APD by preventing the injection of diffused electrons (when 
p>n) and holes (when n>p) into the avalanche region. These issues will be discussed in 
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more details in the following chapter. By increasing the applied reverse bias across the 
structure, the “inactive” stages are slowly depleted of carriers, but the non-symmetry 
characteristic of the field profile is nevertheless preserved as shown in Figure 3-13 for the 
same structure at -20 V. In the case where p > n, the field magnitude gradually increases 
from left to right. The reverse is true in the case where n > p. As seen in Figure 3-13, the 
magnitude of the field throughout the avalanche region varies by almost a factor of four at 
reverse bias of 20V. As a result, certain regions of the structure would reach breakdown 
field sooner than others, and the device’s breakdown voltage is substantially lowered. 
Once the device reaches breakdown voltage, the structure becomes increasingly 
conductive. If the doping imbalance is too high, some regions of the junction will never 
become depleted. Further depletion of these regions will then be limited by the avalanching 
mechanisms leading to breakdown. However, if the doping imbalance is small, the device 
will deplete fully with increasing bias. Figure 3-14 shows the band diagram of such a 
structure where the MQW region is shown to be fully depleted at -20 V in the case where 
p=1.2n. 
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Figure 3-12: Conduction/valence band plot for an unbiased, doped MQW APD where p=1.2n 
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CHAPTER m 


THEORETICAL RESULTS 
3.1 Introduction 

In this chapter, we will present some of the theoretical results that were obtained 
using the Atlas 2-D device simulation framework. Different models were developed for 
the various APD structures using the parameters given in Appendix A for the 
GaAs/AlGaAs material systems. Some examples of the algorithms that were used for the 
simulations are listed in Appendix B. In order to simulate devices with non-uniform band 
structures (i.e., MQWs), the drift-diffusion model with position dependent band structure 
was used. Unlike the hydrodynamic model 26 , the drift-diffusion model neglects “non- 
local” transport effects such as velocity overshoot and energy-dependent impact 
ionization. Velocity overshoot occurs when carriers enter a high-field region where the 
field magnitude exceeds its threshold value. This will cause the carriers to accelerate to a 
higher velocity before relaxing to their equilibrium transport condition. Ignoring these 
effects may have significant impact on submicron structures, but the consequences are 
minor for large devices. It is important to note, however, that the purpose behind our 
modeling was mainly to understand the effect of variations in structure parameters on 
device performance and not to determine exact values of the output variables. In all 
models, Newton's two-carrier method was used for solving Poisson's and the continuity 
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equations. Impact ionization was modeled according to Selberherr(ref), and CV solutions 
were obtained through small signal ac analysis at 1 MHz and with a small signal bias of 
0.03 V. Light IV characteristics were modeled using a 1 mW/cm 2 monochromatic light 
source operating at 632.8 nm, and spectral response simulations were performed with a 1 
W/cm 2 broadband light source ranging from 200 to over 900 nm. 

In order to simplify our models, reduce the number of degrees of freedom, and 
decrease program execution times, the following assumptions were made regarding the 
simulated structures: 

• All structures were assumed to have rectangular geometries having a cross sectional 
area of 75 pm 2 . Therefore, a plane junction approximation was used, and a 1-D 
cartesian coordinate solution to Poisson’s equations was applied. 

• Only SRH and Auger recombination mechanisms were considered. Optical and 
surface recombinations were neglected. 

• The presence of defect/trap centers in bulk materials and at interfaces was neglected. 

• The top (p+) and bottom (n+) regions are uniformly and equally doped. 

• The p and n contacts are perfect ohmic contacts. 

• Doping imbalance in the MQWs is constant throughout an entire structure. 

• Bandgap narrowing effects in AlGaAs are similar to those in GaAs. 


36 



3.2 Simulation Results for a PIN Structure 


We first tested the accuracy of our model by examining the results for the internal 
and external properties of a simple 3 pm thick GaAs PIN structure where the top (p+) and 
bottom (n+) layers were uniformly doped at lxlO 18 cm' 3 . 


3.2. 1 Internal Physical Properties 

Figure 3-1 shows the energy band diagram at zero bias where the energy gap at 
300 K is about 1.43 eV. The built-in voltage is shown to be equal to approximately 1.375 
V. This can be easily verified using equation [2-45] with N A = No = 1x10 cm' , 
ni=1.79xl0 6 cm' 3 , and kT/q=0. 02586 V. The effects of bandgap narrowing due to heavy 
doping is also shown in Figure 3-1. These effects were included in the model as spatial 
variations in the intrinsic carrier concentration 27 : 
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This spatial dependence of n; e results in an adjustment to the electric field profile as 
shown by the presence of the two interface peaks in Figure 3-2. 
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3.2.2 External Properties 


Figure 3-3 shows the simulated dark and light IV characteristics for a PIN APD 
with a 1 pm intrinsic GaAs region. The distributed contact resistance for both the p and n 
contacts was set to lxl 0 4 Q.cm 2 . The device shows a gain of about 4 at a reverse bias of 
38 V. 

Figure 3-4 shows the simulation results for the spectral output of the same PIN 
structure under constant power illumination compared to the output that would be 
expected if no recombination processes were present. The spectral response drops sharply 
around 0.9 pm due to the absorption properties of the material. This long-wavelength 
cutoff is determined by the GaAs energy gap which is about 1.43 eV (~ 867 nm) at room 
temperature. The spectral response data can be used to determine the external quantum 
efficiency of the structure as will be shown in the following chapter. 

Figure 3-5 shows the photogeneration rate in a 3 pm GaAs PIN APD under 400, 
632, & 850 nm illumination. At a wavelength of 400 nm, most of the photogenerated 
carriers are created near the top surface within the first 0.2 pm. At the HeNe wavelength 
of 632 nm, most of the incident light is absorbed in the p + GaAs layer before reaching the 
avalanche region. As the wavelength increases toward the near-infrared part of the 
spectrum, carrier photogeneration in the depletion region becomes increasingly 
important, and the light absorption capacity of the structure is reduced as was 
demonstrated by the spectral output of Figure 3-4. 
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[26] Silvaco International, ATLAS User’s Manual, Device Simulation Software, version 
4.0, June 1995. 

[27] J.W. Slotboom, “The PN Product in Silicon”, Solid State Electronics 20, pp. 279- 
283. 1977. 
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CHAPTER IV 


EXPERIMENTAL APPARATUS 
4.1 Description of Structures 

All of the APD structures were grown by molecular beam epitaxy (MBE) in a 
Varian Gen II system and were fabricated using standard photolithographic techniques. 
Figure 4-1 shows a cross section of an experimental MQW APD where contact rings were 
used to provide a uniform electric field across the entire device. The MQW region in all 
tested devices consisted of 10 to 30 sets of alternating layers of GaAs (500 A) and 
^ 0 . 42 ^* 0 . 58 ^ (500A) with 1000 A periods. All APDs were composed of a 1 mm Be- 

doped (3xl0 18 cm -3 ) p + top layer, and a 1.5 mm Si-doped (3xl0 18 cm' 3 ) n + back layer. 
In the doped-well MQW devices, high electric fields were achieved in the narrow bandgap 
GaAs wells of the avalanche region through the introduction of thin (50-150 A) and highly 
doped (0.5xl0^ 8 -1.5xl0^ 8 cm' 3 ) p + and n + layers. The doped-barrier MQW structures 
were similar with the exception that the doping layers were in the AlGaAs barriers. In the 
undoped MQW structures, no doping was incorporated in the avalanche region, and in the 
conventional APD design (PIN), the MQW region was replaced by a 1 to 2.5 mm intrinsic 
GaAs layer. Pure electron injection was achieved by focusing a 632.8 nm HeNe laser 
inside the circular p-contact ring using a microscope objective. 
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Figure 4-1: Cross Section of a MQW APD 
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4.2 Summary of Capabilities of the APD Characterization Lab 


The purpose of this lab is to measure the various optical and electrical properties 
that determine the output performance of the avalanche photodiode (APD) structures that 
are grown and fabricated at the Quantum Microstructures Lab at GTRI. Various computer 
automated experiments have been developed to study and optimize the different 
parameters of the APDs in order to determine their gain, noise level, breakdown voltage, 
spectral response, and speed. A brief description of some of those experiments and the 
information they provide about the photodiode is presented in the table below. 


Table 4-1: Summary of APD characterization capabilities 


Experiment 

Description 

Information Provided 

Spectral Response 

Responsivity as a function of 
wavelength and gain 

Sensitivity of photodiode 
throughout the optical spectrum 

I-V response 

Diode current output as a function 
of applied bias (10 K - 373 K) 

Photodiode gain and variations 
with incident light intensities 

C-V response 

Diode capacitance as a function of 
applied bias (10 K- 373 K) 

Profile of the carrier concentrations 
vs depletion width in doped APD's 

Noise measurement 

Noise level as a function of laser 
power density (10 K - 373 K) 

Excess noise factor at various 
diode gains and temperatures 

Transient Response 

APD response to a high power 50 
ps laser pulse 

APD bandwidth limitations and 
heterojunction interface quality 


In the following few sections, brief descriptions and illustrative diagrams will be provided 
to explain the operating principles behind each one of the above experiments. 
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4.3 Experimental Techniques 


4.3.1 Spectral Response Measurements 

The spectral response measurements were conducted using a Vi m Jarrell-Ash 
spectrometer and a broadband light source placed near the entrance slit. The light 
component at the exit of the spectrometer was then focused on the top p-regjon of the 
photodiode. The output current of the APD was then measured as a function of 
wavelength between 300 and 900 nm and adjusted for the variations in the lamp’s optical 
spectrum supplied by the manufacturer. By measuring the incident light power and 
neglecting surface reflection, the spectral data can be used to get an estimate of the 
external quantum efficiency of the device. 

4.3.2 Current as a Function of Reverse Bias (IV) 

IV measurements are conducted using a computer-controlled Keithley Source- 

Measure Unit (SMU) which provides the applied voltage and reads the corresponding 
output current of the photodiode. Both dark and light IV measurements are conducted 
and the data is then used to calculate the multiplication gain of the device as a function of 
applied bias. Either electron- or hole-injections can be achieved by focusing the HeNe 
laser beam inside the p-contact ring at the top, or the bottom n-layer as shown in Figure 4- 
1. The operating temperature can be varied between 80 to 400 K by mounting the 
structures inside a Joule-Thomson cooling system supplied by MMR Technologies. A 
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schematic of the optical set-up used to mount and test the APDs is shown in Figure 4-2. 
The vacuum chamber used for low temperature measurements is shown in Figure 4-3. 



Figure 4-2: Experimental Set-up for IV, CV & Noise Measurements 



Figure 4-3 : Vacuum Chamber used for Low Temperature Measurements 
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4.3.3 Capacitance as a Function of Reverse Bias (CV) 


CV Measurements were made using a computer-controlled HP4277A LCZ meter 
which can be interfaced with the same experimental set-up shown in Figure 4-2 and Figure 
4-3. The CV data is then used to calculate the doping concentrations and the depletion 
width profiles of the structures as was explained in Chapter II. 

4.3 .4 Noise Measurements 

Excess noise measurements are made using an HP8568B Spectrum Analyzer 
controlled through an IEEE interface. The APD noise output is measured for a constant 
value of the gain or reverse bias, and is normalized to the noise value corresponding to a 
unity gain. The calculated excess noise factor, F(M), is defined as the “excess” output 
noise resulting from impact ionization processes. It can be measured using either electron 
or hole injection, F e (M) or F P (M). The excess noise factor is then plotted as a function of 
structure gain and the data is compared to McIntyre theoretical curves to determine the 
electron/hole impact ionization ratio. 

4.3.5 Transient Response Measurements 

Time response measurements were conducted using a high-power Hamamatsu 
laser pulser (PLP-03) operating at ~820 nm with a pulsewidth around 50 picoseconds. 
The APD output pulse is detected using a Tekll801 digital oscilloscope w/ a 50 MHz 
SD-32 sampling head. The corresponding set-up is shown in Figure 4-4. Transient 
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response measurements are used to determine transit and diffusion times as well as the 
effect of carrier trappings on overall device bandwidth. 
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Figure 4-4: Experimental Set-Up for Lifetime Response Measurements of APDs 
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4.4 Modeling Capabilities 


Silvaco’s device simulation software uses powerful numerical techniques to solve 
for the various microscopic and macroscopic properties of heterojunction devices. The 
following is an overview of the general capabilities of ATLAS which includes the 
following tools and extensions: 


• ATLAS: 


Supplies general capabilities that are accessible by all device 
simulation products. 


• S-PISCES: 


Simulates silicon devices. 


• BLAZE: Simulates devices fabricated using arbitrary semiconductors 

(including II- VI, III-V, and IV-IV materials), and 
heterojunction devices. 

• GIGA: Adds the ability to perform nonisothermal calculations that 

include the effects of lattice heating and heat sinks. 

• TFT: Allows the simulation of polycrystalline- and amorphous- 

based devices 


• LUMINOUS: Provides capabilities to model optoelectronic devices, 

including sophisticated ray-tracing 

• LASER: Allows the simulation of heterostructure lasers by self- 

consistent solution of the Helmholtz equation for the optical 
field. 


• MIXEDMODE: Offers circuit simulation capabilities that employ numerical 

physically-based devices as well as compact analytical models 

• DEVICE3D: Provides capabilities for three-dimensional device simulation 

• INTERCONNECT3D: Provides capabilities for three-dimensional parasitic 

extraction. 
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• THERMAL3D: Provides capabilities for three-dimensional thermal analysis. 

Using the above tools and packages, one can design programs where the material 
parameters and device structure are defined in the input deck. ATLAS is then used to 
provide a comprehensive set of physical models including -but not limited to- the 
following: 

• DC, AC small-signal and full time dependency analysis. 

• Drift-diffusion transport models. 

• Advanced mobility models. 

• Graded and abrupt heterojunctions. 

• Ohmic, Schottky, and insulating contacts. 

• SRH, radiative, Auger, and surface recombination. 

• Local and non-local impact ionization. 

• Band-to-band and Fowler-Nordheim tunneling. 

• Optoelectronic interactions with general ray tracing. 

• General electronic circuit environment. 

ATLAS supports a large number of semiconductor materials and compounds such as Si, 
GaAs, AlGaAs, ZnS, as well as various metals and insulators such as aluminum, gold, 
SiN, Si02, vacuum, and air. 
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CHAPTER V 


EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 

5.1 Gain/Noise Properties of Doped MOW Junctions 

5.1.1 Introduction 

In this section, a detailed characterization has been made of the external 

properties of both the PIN and the MQW structures. Comparison of the gain properties at 

low voltages between the MQW and conventional APDs showed a direct experimental 

confirmation of a structure-induced carrier multiplication due to interband impact 

ionization. Similar studies of the bias dependence of the excess noise characteristics show 

that the low-voltage gain is primarily due to electron ionization in the MQW-APDs, and 

to both electron and hole ionization in the conventional APDs. For the doped MQW 

APDs, the average gain per stage was calculated by comparing gain data with depletion 

width and carrier profile measurements, and was found to vary from 1.03 at low bias to 

1.09 near avalanche breakdown. These results are in good agreement with theoretical 

28 

models developed by Brennan for similar derivatives of the doped MQW APD. 

5.1.2 Gain Enhancement 

As was previously mentioned, superlattice multiplication APDs are designed 
to outperform bulk multiplication APDs by artificially enhanced ionization through the 
introduction of multiple quantum well layers. This behavior is attributed to the large 
difference in the conduction and valence-band edge discontinuities at the AlGaAs/GaAs 
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interface. To better understand these characteristics, consider once again the energy band 
diagram shown in Figure 1-2. When a “hot electron” enters from the AlGaAs barrier 
layer into a GaAs well, it abruptly gains an energy equal to the conduction band gap 
discontinuity, AEc. The effect is that the electron “sees” an ionization energy reduced by 
AEe with respect to the threshold energy in bulk GaAs (Eth=2.0 eV ) 29 Since the impact 
ionization rate a increases exponentially with decreasing Eth, a large increase in the 
effective a compared to that in bulk GaAs is expected. When the electron enters the next 
barrier layer, the threshold energy in the AlGaAs material is increased by AE C therefore 
decreasing the value of a in the AlGaAs. However, since Og»a» » 0£ajg»a*, the 
exponential dependence on the threshold energy results in an increase in the overall 
average a given by: 

Cl»vg = (cXg»Ai ^ CtA)G»AiV(LG»Ai ^ La1G»Aj) [5*1] 

where L represents the layer thicknesses. 

In contrast, the ionization rate for holes, 3, is not increased substantially due to 
the smaller valence-band discontinuity. This results in a net enhancement in the a/3 
ratio. 30 

The APD devices were characterized under both light and dark conditions 
using current-voltage (I-V), capacitance-voltage (C-V), and noise measurements. Gain 
curves were calculated from the reverse bias I-V measurements performed as a function 
of photon flux. Carrier concentrations and depletion width profiles were determined 
from the C-V data using a one-sided junction approximation. 
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The I-V measurements were taken in the dark and under HeNe laser 
illumination using the experimental apparatus described in section 4.3.2. In order to 
ensure pure electron injection, the laser beam was focused through a microscope 
objective at the center of the 75 jam diameter opening in the upper p-contact ring. 
Breakdown voltages, Vg, were measured to be about 27 V for the doped MQW APD and 
63 V for the conventional APD with corresponding dark currents, measured at 20% of the 
breakdown voltages, of about 10 pA and 100 pA, respectively. The dark current IV plots 
are shown in Figure 5-1 where the low breakdown voltage characteristics of the doped 
MQW APD is demonstrated. This is a result of the high doping present in the junction 
which helps increase the electric field magnitude closer to the its critical avalanche value. 

The C-V measurements were performed at 1 MHz using the apparatus 
described in section 4.3.3. The C-V data (shown in Figure 5-2) was then analyzed to 
calculate the depletion widths and carrier profiles for the two structures The net carrier 
concentration for the conventional APD is shown in Figure 5-3 as a function of the 
calculated depletion width. As the reverse bias is increased, the capacitance decreases to 
0.8 pf, while the depletion width increases to about 2 6 pm prior to breakdown around 
63V. Note that the carrier concentration increases sharply as the depletion edge is 
extended into the doped contact region. Figure 5-4 shows the corresponding plot for the 
doped-well MQW device which had a breakdown voltage of about 27 Volts . This plot 
clearly shows evidence of the depletion of all 10 stages in the MQW APD. The presence 
of the peaks in the carrier profile data is due to the unequal p and n doping concentrations 
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Reverse Bias (V) 


Figure 5-1 : Dark IV plots for the PIN and the doped-well MQW APD 
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Figure 5-3: Carrier concentration profile vs depletion width for the PIN APD. 
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Figure 5-4: Carrier concentration profile vs depletion width for the doped MQW APD. 
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in the wells. This results in partial depletion of the junction as will be shown shortly. 
With applied bias, the depletion width increases, and additional quantum-well stages 
become depleted. This gives rise to ripples in the CV profile resulting from peaks in the 
carrier concentration. Note that the peak positions do not agree with the 1000 A period of 
the MQW structure. This discrepancy is due to the fact that carrier concentration profiles 
were calculated assuming a one sided depletion. 31 This is generally not the case in such 
structures unless there is a large doping imbalance in the junction preventing it from 
depleting both ways. In addition, the spatial resolution of the C-V measurements was 
limited by the Debye length given by 32 , 

l d = fiifoJfF) [5-2] 

which is about 40 A at room temperature for a doping level of n=1.5xl0^ cm'3. The 
Debye length is the distance over which the Coulomb (electrostatic) forces between 
charged layers are essentially screened out. Since the thickness of the doped layers in the 
wells was of the same order of magnitude (50 A), abrupt changes in the doping 
concentration could not be accurately measured. 

The gain curves, calculated from the I-V data, are shown in Figure 5-5 where 
the bias values were normalized by the breakdown voltage of each device to enable 
comparison. Figure 5-5(a) clearly shows the presence of gain in the doped MQW device 
in the low voltage region while the conventional p-i-n structure (Figure 5-5(b)) does not 
show any gain in this regime. This is an indication of a structure-induced carrier 
multiplication resulting from the band discontinuity and the doping in the MQW APD. In 
order to calculate the gain per period in the doped MQW, the carrier profile plot (Figure 
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5-4) was superimposed on the gain curve (Figure 5-5(a)) and the gain was estimated at 
each consecutive carrier concentration minimum as shown in Figure 5-6. The 
corresponding gain values per stage were found to increase from 1 .03 at low bias (one 
depleted stage), to about 1 .09 near breakdown (ten depleted stages). These results are in 
good agreement with theoretical predictions provided by Brennan 28 for similar 
derivatives of the doped MQW APD. 
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Figure 5-5: Gain vs. the ratio of reverse bias to breakdown voltage for the (a)MQW APD 
and (b)PIN APD. 
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Figure 5-6. Plot used for determining gain per stage for a doped MQW APD 
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5.1.3 Dark Current Reduction 


Dark current is one of the main parameters of concern in photodetectors. Dark 
current analysis and the reduction of dark current are very important for high sensitivity 
and low noise applications. Defect centers, heterojunction interface traps, as well as mesa 
surface leakage can generate high levels of excess dark current and reduce minority- 
carrier lifetime. 33 

In a typical PN junction, the overall dark current is the sum of the bulk and 
surface components. The bulk component is usually made up of diffusion, generation- 
recombination, and tunneling currents. The surface component consists of generation- 
recombination, and leakage shunt currents usually formed at semiconductor and dielectric 
interfaces. 34 In a device structure with top p and n contacts such as the APDs used in our 
experiments, there are additional sources of dark current components These are due to 
defect centers at the GaAs/AlGaAs interface and most importantly to surface leakage 
currents along the mesa edge which can contribute significantly to the dark current. 

In this section, it will be shown how substantial the surface leakage component 
can be and how certain growth, processing and surface treatment techniques can be used 
to dramatically lower surface leakage currents by several orders of magnitude. The 
devices that were measured were volume- and delta-doped MQW. In the volume-doped 
MQW structures, the GaAs wells were doped with 50 A (3.0xl0 18 cm'3) adjacent p + and 
n + layers. In the delta-doped APDs, p + and n + layers with a sheet charge density of 1-5 
xlO 12 cm' 2 were introduced separated by undoped spacer layers ranging from 50 to 150 
A. Through careful dopant calibration, the devices were grown such as to achieve full 
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depletion at low bias. After processing the devices into mesa diodes, various surface 
passivation treatments were investigated. These include both plasma ashing in an O 2 
plasma and ammonium sulfide treatments. Through the application of such treatments, a 
decrease in the reverse bias dark current by as much as a factor of 1000 was achieved in 
the low bias region. This can be seen in Figure 5-7 where the dark current is plotted both 
before and after surface treatment by ammonium sulfide. The dark current approximately 
follows a square-root behavior at low to medium reverse bias while at high biases, 
avalanche currents dominate. The rapid increase in the dark current at low reverse bias 
and the large drop in its value achieved by surface treatment are indicative of surface 
leakage. The leakage in these heterojunction mesa diodes was dominated by 
generation/recombination current near the intersection of the mesa surface with the 
GaAs/AlGaAs depletion region. 

As a result of surface treatment, dark currents as low as 1 pA were obtained under 
zero applied bias. In some APDs, the dark currents increased to only 12 pA at 20% of 
breakdown. In addition, these devices exhibited extremely high gains which exceeded 
10,000 in some cases. In most traditional APDs, the presence of high dark currents 
usually presents a limiting factor preventing the further increase in a device’s 
photocurrent gain beyond avalanche breakdown. By reducing the dark currents in these 
devices, it was possible to maintain it at levels well below that of the photocurrent. This 
made it possible to achieve and sustain high levels of gains well beyond breakdown. 
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5.1.4 Excess Noise Reduction 


There are many different types of noise that may be present in an electrical 
device. A brief discussion of the various types of noise mechanisms is presented below. 

5. 1.4.1 Johnson Noise 

Johnson noise is caused by the random motion of thermally energetic electrons in 
resistive materials. Its instantaneous amplitude is not predictable, but the probability of 
its amplitude being within an interval of dV volts is equal to p(V)dV where p(V) is 
expressed by the familiar Gaussian probability function: 

p(y)= Wr‘ v ' n ’' [5 - 31 

where the parameter o is the rms value of the fluctuations and the quantity universally 
accepted to describe the noise output from a resistor, o is bandwidth dependent and is 
expressed as follows: 

ct = ( AkTR,B) m (volts) [5-4] 

where k is Boltzmann’s constant, T is the resistor temperature in K, R« is the resistance in 
ohms, and B is the noise bandwidth in hertz. Johnson noise is “white noise”, that is the 
rms value per unit bandwidth (rms density) is constant from DC to frequencies extending 
into the infrared region. 
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S. 1.4.2 Shot Noise 


Shot noise is the result of random current fluctuations in vacuum tubes and 
semiconductor junctions. It is caused by the random arrival of discrete electron charges at 
anodes, collectors, and drains. The rms value of shot noise is given by: 

/*, = (2 el^BY 2 (A) [5-5] 

where e is the electron charge, I* is the average DC current through the diode, and B is 
the noise bandwidth in hertz. 

5. 1.4.3 Flicker Noise 

Flicker noise is characterized by its spectral composition and for most electronic 
devices, it dominates thermal and shot noise from DC to about 100 Hz. Although flicker 
noise can be detected in virtually all conducting materials with applied power, it seems to 
be most prominent where electron conduction occurs in granular or semiconductor 
devices. For most semiconductor devices, flicker noise is due to surface effects resulting 
in random carrier recombinations at interface traps. Flicker noise exhibits a 1/f" power 
spectrum, with n typically ranging from 0.9 to 1.35. 

5. 1 .4.4 Total Non-multiplication Noise 

Because all the noise sources are considered to be random and uncorrelated, the 
noise power in a system is additive, and the total rms noise is the square root of the sum 
of the squares of each of the three noise sources previously described. The total non- 
multiplication noise output voltage is given by : 

= [4kTR3 + (WR.) 2+ er 2 ] 1/2 volts rms [5-6] 
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5. 1.4.5 Excess Noise 


The excess noise factor is the component of total noise corresponding to 
fluctuations in the process of carrier multiplication in an avalanche photodiode. It is 
defined as the ratio of multiplication-related noise to that of the non-multiplication noise 
defined by equation [5-6]. Excess noise measurements are usually conducted at high 
frequencies where the shot noise is the most dominant non-multiplication term. 
Therefore, all other noise processes are commonly ignored in excess noice computations. 
McIntyre has shown 35 that the statistical nature of the multiplication process adds an 
additional component to the noise which can be included with the shot noise of the APD 
as an excess noise factor. The excess noise factor in the case of pure electron injection is 
given by: 

[5 - 71 

where M is the multiplication factor, and k is the effective electron to hole ionization 
ratio of the APD. The root mean square noise current <i„ 2> can be expressed as: 

<i„ 2 > = 2eIpoM 2 FB [5-8] 

where Ipo is the primary multiplied photocurrent. In other words, the actual photocurrent 
is given as: 

Iph = Ipo * M for Iph » I D (dark current) [5-9] 

A plot of F(M) vs. M from McIntyre’s theory is shown in Figure 5-8 for k’ (=l/k=p/a) 
ranging from 0.001 to 1000. The plots are approximately symmetric on a log-log scale 
about the axis F(M)=M for k and 1/k. At any given gain, lower excess noise is obtained if 
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the carrier with the higher ionization coefficient is injected into the multiplication region. 
For low k’, if the correct carrier is injected, the excess noise can be quite low, with a 
limiting value of 2 for k’=0 at high gain. However, if the wrong carrier is injected, the 
excess noise becomes very high, with the penalty becoming increasingly more severe as 
the disparity between ionization coefficients decreases. Thus, it is important to inject the 
carrier with the higher ionization coefficient into the multiplication region. The lower the 
k’ (or higher the k), the higher the relative difference between a and 3, and the lower the 
excess noise. It is important to note that the McIntyre model is not well suited for 
describing the noise characteristics of MQW devices since it was intended mainly for 
conventional APDs. Better models have been developed by Teich et al. 36,37 and are 
described in the literature. In addition, Marsland ’ and Hayat have recently 
considered the “dead space” between ionization events in their excess noise calculations. 
They concluded that McIntyre’s calculations overestimate the excess noise factor for a 
given k. McIntyre curves were used in our excess noise factor plot for comparison 
purposes in order to clearly illustrate the difference between the noise properties of 
conventional and MQW APDs. Multiplication noise measurements were conducted on 
both APDs using an HP8568B spectrum analyzer set at a 200 kHz center frequency with 
a 10 kHz resolution bandwidth. 
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Figure 5-8: Excess noise factor F(M) vs. M for constant k\ from McIntyre’s equation. 

Excess noise factor measurements were made with a HeNe laser using 
electron injection into the top p + layer. The experimental excess noise factor data is 
shown in Figure 5-9 where the dashed lines represent McIntyre's calculated theoretical 
curves. Figure 5-9(a) for the doped MQW APD clearly shows that for low gains (M<4), 
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the ionization ratio is greatly enhanced (k= 10-50) as compared to that in bulk GaAs 
(k=1.67). This fact is clear evidence of the validity of our previous results for the gain 
values per stage which assume single carrier multiplication at low voltages. At higher 
voltages, however, the value of k is reduced since the holes gain more energy from the 
applied electric field and are more likely to impact ionize 41 . The noise data for the 
conventional APD displayed in Figure 5-9(b) shows the high noise (k~l) characteristics 
of the conventional APD even at low bias voltages. Excess noise factors at higher gain 
values were difficult to obtain since the dark current becomes large at high bias. Note 
that in Figure 5-9(a), the point where the excess noise data break away from the high k 
McIntyre curves corresponds to the breakdown voltage of the doped well APD. In 
addition, at high gains the k ratio for the doped well APD approaches the bulk GaAs 
value of 1.67. This is expected at high fields since the band bending resulting from the 
MQW structure becomes insignificant compared to that induced by the externally applied 
field. 
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Gain 


Figure 5-9: Excess noise factors for both (a) the doped MQW APD, and (b) the PIN APD. 
Dashed lines correspond to McIntyre theoretical curves for k=l,2,5,10, and 50. 
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5.2 Spectral Response Properties 


5.2.1 APD Quantum Efficiency 

The external quantum efficiency of a photodiode is defined as the number of 
electron-hole pairs generated at the output photocurrent per incident photon : 

*L = (V?) / ('V /,v ') [5-7] 

where I p is the photogenerated current as a result of the absorption of incident optical 
power Pop, at a given wavelength. Another related quantity is the responsivity which is 
defined as the ratio of the output photocurrent to the incident optical power: 

<r = A. = A/w [5 . 8] 

P opt hv 1.24 

The quantum efficiency of a photodetector is primarily determined by the absorption 
coefficient a of the material. Figure 5-8 shows the measured intrinsic absorption 
coefficient for several materials used in photodetectors. 1 From that figure, we can see how 
the room temperature absorption for GaAs material drops sharply around 0.9 pm. This 
long-wavelength cutoff wavelength is determined by the GaAs energy gap which is about 
1.43 eV (~ 867 nm) at room temperature. At short wavelengths, the values of a become 
very large, and the radiation gets absorbed very rapidly near the surface where the 
recombination time is short. This will cause the photocarriers to recombine before they are 
collected by the junction region in a photodiode. 
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Figure 5-8 : Optical absorption coefficients for different materials at 77 K and 300 K. 2 
The internal quantum efficiency of a photodiode r^, depends on the wavelength of the 
light as well as the thickness and doping of the absorption material. The absorption 
follows Beer’s law and the internal quantum efficiency can be expressed as = 1- 
exp(ax,b.), where a is the wavelength dependent absorption coefficient as shown in Figure 
5-8 , and x,* is the thickness of the absorbing material. The “absorption length”, 1„ is 
defined as 1/a and gives the amount of material needed so that 1/e of the light would be 
transmitted in the absence of reflections. The external quantum efficiency, , given by 
equation [5-7] , includes the effect of reflection, and can be expressed as (1 - R) Tin*, 
Where R is the wavelength-dependent reflectance of the photodetector. 
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5.2.2 Spectral Response Data 


The long-wavelength behavior of the GaAs material was clearly demonstrated by 
spectral response measurements conducted on a doped PIN APD using the apparatus 
described in section 4.3.1. The spectral response output is shown in Figure 5-7. The 
experimental curve shown in Figure 5-7 was in good agreement with calculated spectral 
response data for similar devices. In order to calculate the experimental quantum 
efficiency, we measured the APD current output using a HeNe laser beam incident inside 
the p + ring with a total power of about 2 pW (inside a circular area with a 75 pm 
diameter). The experimentally calculated quantum efficiency at 632.8 nm was found to be 
about 19% for a doping level of lxlO 18 cm' 3 . Table 5-1 shows the experimental external 
quantum efficiency as a function of the doping in the p+ layer. 


N a in cm' 3 

TW in % 

l.Ox 1018 

18.9 

2.0x 1018 

9.0 

3.5 x 1018 

7.7 


Table 5-1 : Measured quantum efficiencies as a function of doping for a GaAs PIN 3 

The simulated quantum efficiency curve (assuming zero reflection and no surface 
recombinations) is shown in Figure 5-8 for a 3 pm PIN photodiode at zero bias. The 
external source power density was maintained at 0.01 W/cm 2 Using a 75 pm APD with 
an active area of 1.6x10“* cm 2 , the total incident power on the top p surface calculates to 
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be about 1 .6 |iW. This is comparable to the HeNe laser power incident on the surface 
which was used in our quantum efficiency experiment. Notice how the theoretical 
quantum efficiency at 633 is about 27% which is slightly higher than our experimental 
value. This is largely due to the loss of light due to surface reflection which was 
unaccounted for in our simulation. Even though reflection is neglected in the model, the 
maximum external quantum efficiency does not reach 100% due to the presence of various 
carrier recombination mechanisms (SRH, Auger, etc.) which were previously described. 

It is possible to increase the quantum efficiency of the device through the introduction of a 
heavily doped p ++ GaAs top layer which will help create a high-field region to enhance the 
diffusion of photogenerated electron-hole pairs toward the depletion region. Figure 5-9 
shows the calculated improvements in quantum efficiency for various doping differences 
between the 0. 1 pm p ++ layer and the 1 pm p + layer. According to the model, it should be 
possible to increase the quantum efficiency by over a factor of 10 for the high energy part 
of the spectrum (0.2-0.4 pm). In addition, the response throughout the visible spectral 
region becomes more uniform as can be seen in Figure 5-9. 

Note that the introduction of the thin layer does not make a significant difference to the 
quantum efficiency when the top absorption region is heavily doped (3xl0 ,8 cm‘ 3 ). This is 
due to increased carrier recombination in the highly doped 1 pm region which tends to 
reduce the number of carriers diffusing toward the depletion region. 
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Figure 5-10: Theoretical quantum efficiency plot (zero reflection) for a 3 pm GaAs PIN 
(1 pm/1 pm/1 pm) APD where p + = n + =3x!0 18 cm' 3 
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Figure 5-11: Theoretical plots showing possible improvement in quantum efficiency for a PIN 
APD through the introduction of a thin (0.02 mm) top p~ layer (all doping concentrations given 

are in cm' 3 ) 
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5.3 Temperature Dependence 
5.3.1 Impact Ionization vs. Tunneling 


In most materials, the mechanisms of junction breakdown can be attributed to both 
impact ionization and tunneling. One way of determining the dominant process is by 
examining the temperature dependence of the junction’s current-voltage or IV 
characteristics. Since the energy bandgap of GaAs decreases with increasing temperature, 
one would expect the breakdown voltage due to tunneling effect to have a negative 
temperature coefficient. In other word, the breakdown voltage would decrease with 
increasing temperature as shown in Figure 5-10. This is because a smaller applied field 
would be needed to reach the same current levels at higher temperatures. 



Figure $-12: IV characteristics of tunneling breakdown 4 

On the other hand, avalanche breakdown has a positive temperature coefficient where the 
breakdown voltage increases with increasing temperature due to the shorter mean free 
path of carriers at higher temperatures.Increasing T increases the phonon vibrations of the 
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lattice, thereby increasing the probability of premature scattering and reducing the 
ionization coeficients for a given E-field. Other sources of leakage current, such as 
generation-recombination and diffusion also tend to increase with increasing T. The result 
of such an effect is shown in Figure 5-11 by the experimental data obtained for a doped 
MQW APD 

According to Tyagi 5 , the breakdown voltage is related to temperature through the 
following linear relationship: 

V b (T) = V B (T 0 )(\ + fl(T- T 0 )) [5-5] 

Where b>0 in junctions where impact ionization dominate. Such linear dependency was 
shown experimentally to be valid by Forrest et al. 6 in the case of p n junctions. However, 
our experimental measurements has shown that for doped MQW junctions, the data can be 
best fit using a third degree polynomial as shown in Figure 5-12. 
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Figure 5-13: Measured gain as a function of temperature for a doped MQW APD 
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5.3.2 Junction Capacitance and Temperature 


As was shown in equation [2-57], the pn junction capacitance is given by: 

Cj-T. 




2Ma)(jw) 

q N a N d 


-11/2 


[5-6] 


which shows that Cj is proportional to (Vy - Va )’ i/2 where Vbi^kT/q)ln(N A ND/ni 2 ) and 
ni 2 =NcN v exp(-Eg/kT). As was stated in equations [2-14] and [2-15], the effective density 
of states, N c and N v , are proportional to T 3/2 . Therefore, n, decreases with temperature as 
shown in Figure 5-13. The result is an overall positive dependency between Vy and 
temperature. Therefore, V»* increases with temperature which causes the capacitance to 
decrease with increasing temperature as shown in the experimental CV data in Figure 5- 
14. As the reverse bias is increased, V A becomes the dominant term in equation [5-6], and 
the variation of capacitance with temperature becomes less significant. This explains the 
convergence toward a limiting capacitance value at high applied voltages. On the other 
hand, when the temperature is lowered, Vw becomes more dominant, and the capacitance 
decreases at a slower rate with the applied bias (see Figure 5-14). 
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Figure 5-15. Intrinsic carrier densities of Ge, Si, and GaAs as a function of temperature 
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5.2 SPECTRAL RESPONSE PROPERTIES . 
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5.2 Spectral Response Properties 
S.2. 1 APD Quantum Efficiency 


The external quantum efficiency of a photodiode is defined as the number of 
electron-hole pairs generated at the output photocurrent per incident photon : 

TJ txi = (I p /q)/(P opl /hv) [5-10] 

where I p is the photogenerated current resulting from the absorption of incident optical 
power Pop, at a given wavelength. Another related quantity is the responsivity which is 
defined as the ratio of the output photocurrent to the incident optical power: 

U*> A/w [5-11] 

P opl hv 1.24 1 1 

The quantum efficiency of a photodetector is primarily determined by the absorption 
coefficient a of the material. Figure 5-11 shows the measured intrinsic absorption 
coefficient for several materials used in photodetectors. 42 From this figure, we can see that 
the room temperature absorption for GaAs material drops sharply around 0.9 pm. This 
long-wavelength cutoff wavelength is determined by the GaAs energy gap which is about 
1.43 eV (~ 867 nm) at room temperature. At short wavelengths, the values of a become 
very large, and the radiation gets absorbed very rapidly near the surface where the 
recombination time is short. This will cause the photocarriers to recombine before they are 
collected by the junction region in a photodiode. 
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Figure 5-10: Optical absorption coefficients for different materials at 77 K and 300 K. 43 
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LIGHT PENETRATION DEPTH 1/a ifim) 




The internal quantum efficiency of a photodiode r|int depends on the wavelength of the 
light as well as the thickness and doping of the absorption material. The absorption 
follows Beer’s law and the internal quantum efficiency can be expressed as rim, = 1- 
exp(otx,b.), where a is the wavelength dependent absorption coefficient as shown in Figure 
5-10, and x,t» is the thickness of the absorbing material. The “absorption length”, L, is 
defined as 1/a and gives the amount of material needed so that 1/e of the light would be 
transmitted in the absence of reflections. The external quantum efficiency, r|ext , given by 
equation [5-10] , includes the effect of reflection as well as the various carrier 
recombination mechanisms. 

5.2.2 Spectral Response Data 

The long-wavelength behavior of the GaAs material was clearly demonstrated by 
spectral response measurements conducted on a doped PIN APD using the apparatus 
described in section 4.3.3. The spectral response output is shown in Figure 5-11. The 
experimental curve shown in Figure 5-1 1 was in good agreement with calculated spectral 
response data for similar devices. In order to calculate the experimental quantum 
efficiency, we measured the APD current output using a HeNe laser beam incident inside 
the p + ring with a total power of about 2 pW (inside a circular area with a 75 pm 
diameter). The experimentally calculated quantum efficiency at 632.8 nm was found to be 
about 19% for a doping level of lxlO 1 * cm' 3 . Table 5-1 shows the experimental external 
quantum efficiency as a function of the doping in the p + layer. 
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Table 5-1: Measured quantum efficiencies (at 633 nm) as a function of doping for a GaAs PIN* 4 


N a in chi ' 3 

Tl„t in % 

foxio 11 

18.9 

mmm 

9.0 

3.5 x 10 18 

7.7 


The simulated quantum efficiency curve (assuming zero reflection and no surface 
recombinations) is shown in Figure 5-12 for a 3 tun PIN photodiode at zero bias. The 
external source power density was maintained at 0.01 W/cm 2 . Using a 75 pm APD with 
an active area of 1.6x1 O' 4 cm 2 , the total incident power on the top p surface is calculated 
to be about 1.6 |iW. This is comparable to the HeNe laser power incident on the surface 
during the quantum efficiency experiment. Notice how the theoretical quantum efficiency 
at 633 nm is about 27% which is considerably higher than our experimental value. This is 
largely due to the loss of light due to surface reflection and to surface recombination 
mechanisms which were unaccounted for in our simulation. Even though reflection is 
neglected in the model, the maximum external quantum efficiency does not reach 100% 
due to the presence of various carrier recombination mechanisms (SRH, Auger, etc.) 
which were previously described. 

It is possible to increase the quantum efficiency of the device through the 
introduction of a heavily doped p ++ GaAs top layer which will help create a high-field 
region to enhance the diffusion of photogenerated electrons toward the depletion region. 
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Figure 5-13 shows the calculated improvements in quantum efficiency for various doping 
differences between the 0. 1 nm p~ layer and the 1 |im p + layer. According to the model, it 
should be possible to increase the quantum efficiency by over a factor of 10 for the high 
energy part of the spectrum (0.2-0. 4 |im). In addition, the response throughout the visible 
spectral region becomes more uniform as can be seen in Figure 5-13. 

Note that the introduction of the thin layer does not make a significant difference 
to the quantum efficiency when the top absorption region is heavily doped (3xl0 18 cm‘ 3 ). 
This is due to increased carrier recombination in the highly doped 1 (im region which 
tends to reduce the number of carriers diffusing toward the depletion region. 


90 




91 







Figure 5-12: Theoretical quantum efficiency plot (zero reflection) for a 3 pm GaAs PIN 
(lpm/ lpm/ lpm) APD where p + = n + =3xl0 18 cm* 3 . 
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Quantum Efficiency (%) 



Figure 5-13: Theoretical plots showing possible improvement in quantum efficiency for a PIN 
APD through variations in the doping concentrations and the introduction of a thin 
(0.02 pm) top p 4-4 layer (all doping concentrations given are in cm' 3 ) 
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5.3 Time Response Characteristics 


The requirements of high bit rate ( > 4000 Mb/s) lightwave communication and 
image processing systems have necessitated the development of fast photodiodes that have 
higher sensitivity than the PIN detector. Improvement in sensitivity while maintaining wide 
bandwidths can only be provided using an amplification mechanism within the 
photodetector itself as is the case in avalanche photodiodes. For optimum operation, an 
APD must meet the following criteria 45 : (1) the electric field in the multiplication region 
must be high enough to produce sufficient gain; (2) the electric field in the absorbing 
region must be low enough so that the tunneling component of the dark current is 
negligible; (3) the depletion region must extend far enough into the absorbing region so 
that diffusion effects are negligibly small. These requirements impose rather severe 
constraints on the doping concentrations and thickness of the epitaxial layers of the device. 
Transient response measurements can provide valuable information on the speed 
performance and the various factors that affect the bandwidth of APDs Since carrier 
diffusion plays a major role in determining the time response of a photodetector, one 
would expect the speed to depend greatly on the depletion characteristics of the APD. The 
more depleted a structure is, the shorter the distance the carriers will have to diffuse, and 
the faster the response time will be. This will be shortly demonstrated with experimental 
time response data. But first, in order to better understand the results that were obtained, 
a brief discussion the various physical effects that limit the frequency response of a 
photodiode will be presented. 
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5.3.1 Bandwidth Limitations of Photodetectors 


Considerable work has been done in the literature on the study of photodetector 
bandwidth limitations [1-4]. To summarize the results, the bandwidth of a Separate 
Absorption and Multiplication Regions (SAM) APD is determined by five physical effects: 

1. Transit Time: This is the time it takes the generated carriers to travel through the 
depleted region under the effect of the electric field. There are two types of transit times in 
an APD. The primary carrier transit time corresponding to the photogenerated carriers, 
and the secondary carrier transit time required for the multiplied carriers of opposite type 
to retrace the steps of the primary ones. Transit times for electrons and holes (x e and Th) 
are usually calculated using the ratio of the distance traveled and the "saturation” velocity 
of the appropriate carrier. 

2. Carrier diffusion time: In the undepleted regions of the device, carrier transport must 
take place by diffusion rather than drift. Because of the absence of electric field in the 
absorption layer of a PIN APD, the photogenerated carriers must diffuse in order to reach 
the avalanche region. This results in the slowing of the device’s response. An oscilloscope 
trace of the transient output of such a device would show both a “fast” and a “slow” 
component. The fast component is due to carrier drift, and the slow one, referred to as the 
“diffusion tail”, is due to diffusion from the undepleted regions of the device. 
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3. RC time constant: There is a fundamental limit on bandwidth due to the capacitive 
transient charging effects which arise from the depletion region capacitance of the device 
and the combined resistance R of the load and the device. 

4. Hole trapping: In heterojunction APDs, there is a possibility of carrier delay caused by 
traps present at the heterojunction interface. This effect is related to the abruptness of the 
heterojunction, the barrier height, the temperature, and the effective mass of the carrier. 
Because the effective mass of holes is larger by an order of magnitude than that of 
electrons, trapping is more likely to occur for holes than electrons. This phenomena is 
known as "hole trapping" and it can be minimized through the use of graded composition 
layers instead of abrupt heterojunctions. 

5. Avalanche buildup time: For single carrier ionization, one only needs to consider the 
transit time through the multiplication layer. For dual carrier ionization, however, there is 
a feedback process that introduces a time delay through the multiplication region. This is 
called the avalanche buildup time. In an APD, there is a buildup time x iV *i associated with 
the avalanche gain process which tends to limit the time response of the photodetector. 
The primary avalanche build-up time for electron initiated multiplication is : 

[5-13] 

V * + V P 0 0 

where K^up is a correction factor; t^, and Up are the electron and hole velocities. The 
physical origin of is the electron/hole displacement current which arises from the 
space-charge induced E-field resulting from the motion of carriers. 46 The closer the value 
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of a is to that of P, the more secondary carriers are generated, and the higher the 
avalanche build-up time as can be seen from equation [5-13], 

5.3.2 Overall Photodetector Bandwidth 

In the case of a PIN photodetector where absorption takes place in the junction, 
the basic limitations to the response time are due to the RC and the transit times of the 
primary carriers. The overall PIN time constant is usually approximated by the square root 
of the sum of squares of the RC and transit time constants: 

T pifl 2 = Max(Th, x e ) 2 + Trc 2 [5-14] 

As was previously mentioned, in an APD, there are two different transit times arising from 
the primary carriers traveling to, and secondary carriers traveling from, the multiplication 
region. In addition, there is the avalanche buildup time (proportional to gain) which is 
proportional to the multiplication process. The actual APD frequency response is a 
complicated function of all of these processes. Hollenhorst 47 and Roy 48 have developed 
complicated transfer functions and matrix expressions to estimate the time constants for 
arbitrary structures. For approximation purposes, the RC time constant is usually treated 
as being non-correlated with the rest of the time constants. In addition, the primary hole 
transit, hole trapping, avalanche buildup and secondary electron transit events can be 
assumed to occur in series, one following the other. In this case, the sum of squares can be 
used to approximate the total time constant as well. 

In the following section, we will present some of the experimental data and attempt 
to provide the proper interpretations as they relate to the processes described above. 
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5.3.3 Experimental Results 


Figure 5-14 shows the pulse response for an unbiased doped MQW device with a 
2.5 ji.m MQW region. The APD was mounted on a 50 GHz Tektronix sampling scope and 
was excited with an 810 nm 50 ps laser pulse. As can be seen from the figure, the unbiased 
APD output pulse has a rise time of about 317 ps, a fall time of 2.5 ns and a full width at 
half max (FWHM) of about 1.4 ns. The oscilloscope trace is a convolution of the 50 ps 
gaussian laser pulse with the output response of the APD. The fast rise time component 
typically follows the relaxation oscillation of the laser pulse. The falling edge of the pulse 
shows a "fast" and a "slow" component. The slow component at the trailing edge is 
usually attributed to either charge trapping at interface states or diffusion of carriers in the 
undepleted regions of the structure. 49 Diffusion will limit the speed of the device as long as 
there are undepleted regions in the structure and a separate absorption layer is being used. 
At high bias, trapping is no longer an issue, and the device response is limited by the 
transit time and the RC time constant. In the following, it will be demonstrated that such a 
slow response is due largely to diffusion effects in the partially depleted APD structure. 

Figure 5-15 shows the response of the above APD under bias (low gain). The fall 
time and the FWHM have now dropped to 819 and 952 ps, respectively, corresponding to 
about a 32% increase in the speed of the device. If the bias is increased further, as shown 
in Figure 5-16, those values drop to 570 and 593 ps, respectivley, with a speed increase of 
about 58%. The large dependence of the pulse’s width and tail on the applied bias is a 
clear indication of a diffusion-limited time response. 
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Figure 5-14: Oscilloscope trace of the pulse response of a doped 2.5 jim MQW APD 
under no bias 
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Figure 5-15: Oscilloscope trace of the pulse response of a doped 2.5 pm MQW APD 
under low bias 
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Figure 5-16: Oscilloscope trace of the pulse response of a doped 2.5 MQW APD near 
breakdown. 

The diffusion-limited effect can be verified further by examining the depletion 
width characteristics of the structure obtained from CV measurements. This is shown in 
Figure 5-17 where the zero bias depletion width of the MQW and PIN structures are 
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about 0. 1 (im and 1 .4 tun respectively. These values can be used to estimate the drift time 
constants from Xdrin = Wdepie»«j/v M t where W deleted is the width of the depleted region, and v ul 
is the saturation velocity taken to be in the order of 10 7 cm/s. The drift time constants for 
the MQW and PIN structures were estimated to be 1 ps and 14 ps, respectively. Similarly, 
the ratio of the diffusion time constants between the MQW and the PIN structures is 
directly related to the ratio of the undepleted widths (T*f^w/Tdifl(piN)= 
W^KfcpfMQw/W uodep^iN) = 2.2). Using that ratio and the sums of squares approximation, the 
diffusion time constants for the MQW and PIN structures were calculated to be 1.5 ns and 
0.7 ns, respectively. This result shows that the MQW structure is largely diffusion-limited 
due to the presence of a large undepleted region. As will be shown in Section 5.5, this 
behavior is largely due to a mismatch in the doping balance between the n and p doping 
layers in the MQW structure. On the other hand, the PIN APD shows a much faster time 
response (Figure 5-18) due to the fact that the structure is largely depleted even at zero 
applied bias (Figure 5-17). Therefore, the time response limitations for this structure are 
mainly due to 1) diffusion time in the top p + and bottom n + layers, and 2) transit time in the 
intrinsic field region of the structure. Diffusion time in the cap layers can be optimized by 
varying the thicknesses as well as the doping concentrations. The transit time can be 
shortened by increasing the field (applied bias) across the junction. This, however, begins 
to creates an additional delay near breakdown due to the increase in the avalanche buildup 
time. 

The "ringing effect" seen in the oscilloscope trace following the output pulse was 
due to the impedance mismatch between the APD circuit and that of the oscilloscope 
sampling head. With applied bias, the impedance of the APD changes due to the increased 
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conductivity of the structure. This increased the impedance mismatch in the circuit 
resulting in reflections of the output pulses as illustrated in Figures 5-15 and 5-16. Note 
that the relative magnitude of the pulses in the case of the MQW APD is not 
representative of the gain of the device since a variable resistor was used in an attempt to 
match circuit resistance and thus limited the voltage applied at the oscilloscope. 
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Figure 5-17: Depletion width plot showing the increase in the depletion region with 
applied bias for a largely undepleted doped (2.5 pm) MQW APD. 
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Figure 5-18: Oscilloscope trace of the pulse response of an undoped 2.5 jim PIN APD with no 
applied bias. 

5.4 Temperature Dependence 

5,4, 1 Impact Ionization vs. Tunneling 
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In most materials, the mechanisms of junction breakdown can be attributed to both 
impact ionization and tunneling. One way of determining the dominant process is by 
examining the temperature dependence of the junction’s IV characteristics. Since the 
energy bandgap of GaAs decreases with increasing temperature, one would expect the 
breakdown voltage due to the tunneling effect to have a negative temperature coefficient 
leading to a decrease in breakdown voltage with increasing temperature as shown in 
Figure 5-19. This is because a smaller applied field would be needed to reach the same 
current levels at higher temperatures. On the other hand, avalanche breakdown has a 
positive temperature coefficient where the breakdown voltage increases with increasing 
temperature due to the shorter mean free path of carriers at higher temperatures. 
Increasing T increases the phonon vibrations of the lattice, thereby increasing the 
probability of premature scattering and reducing the ionization coefficients for a given E- 
field. Other sources of leakage current, such as generation-recombination and diffusion 
also tend to increase with increasing T. The result of such an effect is shown in Figure 5- 
20 by the experimental gain data obtained from the IV curves of a doped MQW APD. 
According to Tyagi 50 , the breakdown voltage is related to temperature through the 
following linear relationship: 

v b (T) = F fl (r 0 >(i +/?(r- r 0 » [5-5] 

where P>0 in junctions where impact ionization dominates. Such a linear dependence was 
shown experimentally to be valid by Forrest et al. M in the case of p + n junctions. However, 
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Figure 5-20: Measured gain as a function of temperature (in K) for a doped MQW APD. 


our experimental measurements have shown that for doped MQW junctions, the data can 
be best fit using a third degree polynomial as shown in Figure 5-21 . 


108 



5.4.2 Junction Capacitance and Temperature 


As was shown in equation [2-57], the PN junction capacitance is given by: 


Cj 


K s £qA 

2K t e 0 (V bi -V A )(N- + N;) ' 

<7 A r A N* D . 


[5-6] 


which shows that Cj is directly proportional to [N a 7v1d + /(Na' + N D + )] 1/2 The number of 
ionized donors and acceptors are given by: 43 


N + D =N D 


1- 



n: 


N A 


l + Sa ex PI 


jjv A jv 
kT 


[5-7] 


[5-8] 


where gp is the ground state degeneracy of the donor impurity level and is equal to 2 since 

a donor level can accept one electron with either spin or can have no electron. On the 
other hand, g A is the ground-state degeneracy factor for acceptor levels and is equal to 4. 

This is because in GaAs as well as in Ge and Si, each acceptor impurity level can accept 
one hole of either spin and the impurity level is doubly degenerate as a result of the two 
degenerate valence bands at k=0. 

Therefore, by examining equations [5-6]-[5-8], it is clear that the capacitance is 
expected to decrease exponentially with decreasing temperature following the decrease in 
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Figure 5-21: Experimental data and fits for a doped-well MQW APD. 
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the ionized donor and acceptor densities. This behavior is illustrated in the experimental 
CV data in Figure 5-22 where the capacitance at low bias decreases exponentially toward 
a limiting value of about 2.8 pF in the case of a doped-well MQW structure. As the 
reverse bias is increased at a given temperature, the depletion width increases causing the 
capacitance to drop toward 2.8 pF corresponding to the capacitance value for maximum 
depletion of the structure. 
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Figure 5-22: Experimental CV curves as a function of temperature for a 
doped-well MQW APD. 
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A COMPREHENSIVE ANALYSIS OF THE PHYSICAL PROPERTIES OF 


ADVANCED GaAs/AlGaAs JUNCTIONS 

Hicham M. Menkara 
Directed by Dr. Christopher J. Summers 

In this work, a detailed experimental investigation and analysis were performed of the 
physical properties of advanced semiconductor junctions. The analysis includes a study of (1) the 
difference in the structure-induced multiplication gain between doped GaAs/AlGaAs MQW and 
PIN junctions, (2) the effect of variations in the doping profiles on the properties of doped MQW 
structures, (3) the effect of surface treatments on the dark current and gain characteristics, (4) the 
spectral and time response limitations of the structures, and (5) a complete modeling of the 
junction physics for the different structure types. All of these investigations were performed by 
conducting experimental measurements and theoretical simulations on new avalanche photodiode 
(APD) structures with built-in intrinsic (PIN), doped and undoped MQW structures. 

In order to properly interpret the experimental data provided by the measurements, one 
needs to understand the correlation between such data and the physical parameters used in 
designing the structures. To accurately determine such a correlation usually requires the ability to 
grow and fabricate a large sample of structures produced under very similar conditions. 
Unfortunately, this is not very practical or even possible during material growth and fabrication. 
Therefore, developing theoretical models which accurately predict the relationships between the 
input and the output parameters is essential to understanding the physics behind the data. 



Most of the experimental data will be presented and analyzed for the first time in the 
GaAs/AlGaAs material system. The experimental results were compared to theoretical models, 
and were used to demonstrate, for the first time, the impact of the doping imbalance throughout 
the structure on the optical and electrical characteristics of a doped MQW structure. These 
models accurately predicted most of the external behavior displayed by these structures during 
experimental testing. In addition, various surface treatment techniques which enabled a dramatic 
reduction in the reverse bias dark current by as much as a factor of 1000 will be discussed. 
Furthermore, a new technique will be presented for improving the quantum efficiencies of these 
structures, and its effectiveness was verified through theoretical models. 



S-5 EFFECT OF VARIATIONS IN THE DOPING PROFILES 
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5.5 Effect of Variations in the Doping Profiles 


The purpose of the following analysis is to use both theoretical and 
experimental evidence to determine the impact of doping imbalance and symmetry on the 
physical and electrical characteristics of doped MQW APDs. Theoretical models have 
been developed to calculate the electric field, valence and conduction band profiles, CV 
profiles, as well as carrier concentration versus depth profiles. Our models showed a 
strong correlation between the p- and n-doping balance inside the GaAs wells and the 
number of depleted stages and breakdown voltage of the APD. A periodic doping 
imbalance in the wells has been shown to result in a gradual increase (or decrease) in the 
electric field profile throughout the device which gave rise to partially depleted devices at 
low bias. The MQW APD structures that were modeled consisted of the standard 
structure with a 1 pm doped-well MQW region. These simulation results showed that in 
an APD with nine doped wells, and where the 50 A p-doped layer is off by 10% 

compared to the n-doped layer (p=1.65xl0 18 cm'3, n=1.5xl0 18 cm'3), half the stages 
were shown to be undepleted at low bias which was a result of a reduction in the E-field 

near the p + cap layer by over 50% from its value in the balanced structure. Experimental 
CV and IV data on similar MBE grown MQW structures have shown very similar 
depletion and breakdown characteristics. The models have enabled a better interpretation 
of the experimental data and relate some of the observed peculiarities in the IV and CV 
curves directly to the doping profile in the MQW structure. 
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5.5.1 Introduction 


As was described in Chapter IV, various characterization techniques have been 
devised to analyze and understand the optical and electrical properties of APDs 53 . These 
include various experiments such as IV measurements which are used to determine the 
gain properties of the device, and CV measurements which are used to calculate carrier 
concentration versus depletion width profiles. However, the data obtained using such 
experimental techniques are not always easy to interpret and relate back to the physical 
processes taking place inside the structures. In addition, the practical limitations inherent 
in the growth and fabrication of large quantities of devices with different structural 
designs add even more complexity to the problem because of the large number of 
variables involved in the process. 

In what follows, a more practical approach will be presented to analyzing the 
experimental data obtained using IV and CV experiments and specifically those relating 
the doping profile characteristics to device properties. Accurate theoretical models of 
MQW APD structures have been developed using AtlasII, Silvaco’s two-dimensional 
device simulation framework. These models were used to provide graphical 
representations of the spatial variations of the electric field across the biased structure, as 
well as conduction and valence band diagrams of the GaAs/AlGaAs MQW structure 
before and after breakdown. In addition, avalanche breakdown simulations and small 
signal ac analysis were used to extract IV and CV curves in order to compare the data 
from the models to those obtained directly from our experimental devices. Both electron- 
and hole-injected photocurrent solutions were obtained by simulating a 632.8 nm 
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monochromatic light source with spot power of about 1 W/cm 2 incident on the devices’ 
front and back surfaces. CV solutions were obtained using small signal analysis at a 
frequency of 1 MHz and with a signal magnitude of 0.03 V. All our analysis were 
conducted using Newton’s two-carrier method 13 and the generation rate of electron-hole 
pairs due to impact ionization was modeled according to Selberherr 54 . 

5.5,2 Theoretical Results 

The APD structure used in this model consisted of a top and bottom p + and n + 
doped (3xl0 18 cm* 3 ) GaAs layers with thicknesses of 1 pm. The middle region was made 
up of 10 periods of alternating layers of GaAs (500 A) and Al 0 42 Gao 58 As (500A). The 
GaAs wells were similarly doped with p-i-n layers whose thicknesses and doping 
concentrations were treated as variable parameters for the purpose of our study. When a 
reverse bias is applied, the combined effect of the applied electric field, the built-in field, 
and the conduction band offset enhances the ionization process of electrons in the GaAs. 
The holes, on the other hand, are subjected to a smaller valence band discontinuity and 
therefore gain less energy than the electrons. 

In this study, the widths of the p and n doping layers were held constant at 50 A, 

and that of the intrinsic layer at 100 A. The doping imbalance ( I p \-\»\ ) was varied 

I "I 

between zero and 100 percent. Figure 5-23 shows the corresponding CV plots obtained 
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Figure 5-23: Comparison of theoretical CV data obtained for the same APD MQW structure 
where the doping mismatch in the wells was varied between 0% and 100%. 
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for similar devices with 0%, 33.3%, and 100% doping imbalance. As expected, the 
capacitance of the device is lowest when the p and n doping are perfectly matched since 
the net carrier concentration throughout the device is reduced to zero. However, the 
avalanche breakdown voltage as depicted by the IV curves in Figure 5-24, seems to be 
highest when p is equal to n. This is due to the fact that a doping mismatch would result 
in a gradual increase of the electric field throughout the device which would cause impact 
ionization to take place at a lower bias point (see E-field profile in Figure 5-25). 
Therefore, a large doping imbalance would actually lower the bias at which breakdown 
occurs. This, however, comes at the expense of a large undepleted region which could 
limit the quantum efficiency and severely hurt the time response characteristics of the 
photodiode. 

Figure 5-26 shows the calculated carrier concentration versus depletion width 
profile for similar APD structures with 0%, 10%, 20%, 33.3%, and 100% doping 
mismatch. In the case where p=n=1.5xl0 18 , it can be seen that the device is folly depleted 
at zero bias. The number of undepleted stages begins to increase when increasing the 
offset between p- and n-doping. In the case where the p-doping is twice that of n, only 
about 20% of the device is depleted at zero bias. In order to better understand the effect 
of the doping imbalance on the MQW structures, it is helpful to examine the valence and 
conduction band diagrams (Figure 5-27), as well as the electric field spatial profile shown 
in Figure 5-25. Devices corresponding to 0%, 33%, and 100% doping imbalances are 
modeled at a reverse bias of 20 V. As is seen from the two figures, the electric field is 
uniformly symmetric, and the MQW region is equally depleted in the case where p=n. 
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However, as p gradually increases, the electric field becomes progressively lower near 
the top p-layer which results in non-uniform depletion of the MQW structure. The effect 
of such non-uniform depletion on the device’s photocurrent can be clearly seen in the IV 
plots shown in Figure 5-24. In the case of electron injection, the photocurrent at zero bias 
is about four orders of magnitude lower than that for hole injection. Gradually, the 
electron-injected photocurrent increases as the device is depleted until it reaches about 
6x1 O' 8 A corresponding to that of the hole-injected photocurrent. Therefore, a doping 
imbalance where p>n can greatly reduce the device’s external quantum efficiency in the 
case of electron injection. Such an effect is not as pronounced in the case where n>p due 
to the smaller valence band discontinuity faced by the injected holes. 
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Figure 5-24: Comparison of theoretical light IV data obtained for the same APD MQW 
structure where the doping mismatch in the wells was varied between 0% 
and 100% 
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Figure 5-26: Calculated doping profiles versus depletion width using theoretical CV data obtained 
from AC analysis on similar MQW APDs with different doping imbalances in wells. 
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Figure 5-27: Theoretical band diagrams of the same MQW APD structures with different doping 
imbalance in the wells (V = -20 V) 
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5.5.3 Experimental Results 


CV measurements were performed on all devices at 1 MHz using the 
previously described set-up. The CV data was then analyzed to calculate the depletion 
widths and carrier profiles for the structures. Figure 5-28 shows the net carrier 
concentration plots for four doped-well MQW devices labeled as APD1 through APD4. 
All devices have similar geometries except for the p- and n-doping in the wells which 
was varied between 0.5xl0 18 and 1.5xl0 18 cm* 3 . It is interesting to see that even though 
APD1 and APD2 were expected to have very similar properties, their CV and carrier 
concentration profiles were quite different. APD1 was almost fully depleted at zero bias, 
while APD2 was not and only reached full depletion right before breakdown. Note that 
the peak positions in the carrier profile of APD2 do not quite agree with the 1000 A 
period in the MQW structure and with the doping profiles obtained using our models. 
This discrepancy is due to the fact that carrier concentration profiles calculated from the 
experimental CV data assumed a one sided depletion 55 which apparently does not hold 
true for the experimental devices. Other sources of error in the experimental data result 
from the inability to accurately account for parasitic capacitance between the devices and 
the metal contacts and bonding wires in the measurement system. In addition, note how 
the average net carrier concentration in the experimental doping profile gradually 
increases up to the top GaAs well where it then drops indicating that the doping 
imbalance is not the same throughout the structure. 
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The spatial resolution of the CV measurements (both experimental and theoretical) was 
limited by the Debye length (equation [5-1]) which is about 40 A at room temperature 

for a doping level of n=1.5xl0 18 cm" 3 . Since the thickness of the doped layers in the 
wells was of the same order of magnitude (50 A), abrupt changes in the doping 
concentration could not be accurately measured. Therefore, it is generally difficult to 
relate the apparent carrier concentration obtained from the experimental devices to the 
actual doping imbalance in the wells. However, using our theoretical carrier profile data 
shown in Figure 5-26 where the actual doping imbalance is fully known, it is possible to 
estimate the actual doping mismatch in every doped layer in the experimental MQW 
device by superimposing both the experimental and theoretical data. By examining Figure 
5-26 and Figure 5-28, we can roughly conclude that the average doping imbalance in 
APD1 is far less than 10% which resulted in full depletion at zero bias, while that in 
APD2 is between 30% and 40% where full depletion of the MQWs was achieved near 
breakdown around 27 V. In the case of APD3 and APD4, the situation was quite 
different. Apparently, the p- and n-doping mismatch was so large (~ 200 %) to the point 
where only partial depletion of two wells was accomplished before avalanche breakdown. 
By examining the electron injected photocurrent curves in Figure 5-29, we can easily 
conclude that for both APD1 and APD2, the doping mismatch is such that n > p, while in 
APD3 and APD4, the situation is reversed. This can be clearly understood by comparing 
the light IV data in Figure 5-29 to the theoretical curves shown in Figure 5-24. The 
gradual increase in the electron injected photocurrent in APD3 and APD4 is an indication 
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of trapping of injected electrons by the AlGaAs barriers near the p-layer where the device 
is undepleted. Therefore, according to the models shown in Figure 5-26, the average p- 
doping in the wells must be larger than that of n. In the case of APD1 and APD2, no 
trapping seems to take place since the low bias photocurrent is much higher and relatively 
flat. Therefore, these two devices have undepleted regions near the n-layer which 
indicates that n > p. 
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CHAPTER VI 


CONCLUSIONS 

Throughout this work, a detailed comparison of the gain and noise characteristics 
of a conventional and a doped well MQW APD was presented. The data obtained 
demonstrated a direct experimental evidence of structure induced preferential 
multiplication of electrons over holes. For the doped MQW APDs, the average gain per 
stage was calculated by comparing gain data with carrier profile measurements, and was 
found to vary from 1.03 at low bias to 1.09 near avalanche breakdown. This is in contrast 
to conventional PIN structures which show no gain in this regime. It was also shown that, 
as the bias was increased, the effect of the structure became less pronounced, and the 
MQW device was reduced to a conventional PIN structure. Similar studies of the bias 
dependence of the excess noise characteristics show that the low-voltage gain is primarily 
due to electron ionization in the MQW APDs, and to both electron and hole ionization in 
the PIN APDs. Our measurements of the doped MQW APD clearly showed that for low 
gains (M < 6), the ionization ratio is greatly enhanced (k = ct/p = 10 - 50) as compared to 
that in bulk GaAs (k = 1.67). At higher voltages, however, the value of k is reduced since 
the holes gain more energy from the applied electric field and are more likely to impact 
ionize. 

It was also observed that surface recombination has a significant on the dark 
current behavior of an APD. The resulting leakage currents can have dramatic 
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consequences on the sensitivities and attainable gain levels in a photodiode. As a result of 
surface treatment, dark currents at low bias were reduced to as low as 1 pA. The result of 
this reduction in dark current was manifested in the structures’ high gain performance 
which exceeded 10,000 in some APDs. By being able to reduce the dark currents, it was 
possible to maintain dark current levels well below those of the photocurrents. This made 
it possible to achieve and sustain high levels of gains well beyond the onset of junction 
breakdown. 

The spectral response and quantum efficiencies for some of the structures were 
also calculated and modeled. The experimental data were very consistent with the 
theoretical models. The quantum efficiencies of the fabricated structures were relatively 
low because of the loss of photogenerated carriers due to recombination mechanisms in 
the diffusion layer. It was shown how it was theoretically possible to significantly 
increase the quantum efficiency of the devices through the introduction of a heavily 
doped p ++ GaAs top layer. Such layers help to create a high-fteld region that will enhance 
the diffusion of photogenerated electrons toward the depletion region. 

An investigation was also made into the impact of doping imbalances in doped- 
well MQW APDs on device IV and CV characteristics and how such an imbalance would 
affect the depletion properties of the APDs. Our theoretical models were in full 
agreement with the observed experimental data and have provided a good understanding 
of the physical processes that take place inside a doped MQW APD. These models have 
been used to interpret experimental IV and CV data and to determine the extent of 
depletion in APD devices. How these parameters were affected by the p- and n-doping 
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imbalance in the structure was also determined. The model predicted that a doping 
mismatch as small as 10% could reduce the depletion layer by as much as 50%. It was 
also shown how a large doping imbalance would cause the device to quickly reach 
avalanche breakdown in the depleted layers and would prevent full depletion. 

The presence of undepleted regions was also shown to be one of the major causes 
of the slow time response in avalanche photodiodes. Partial depletion gave rise to a 
diffusion-limited transient response in doped MQW structures. This was demonstrated 
experimentally to be the case by examining the change in diffusion tail of the output 
pulse response of the devices as a function of applied bias. Fully depleted PIN structures 
showed a fast time response even at zero applied bias. The relationships between the 
depleted (undepleted) widths and the drift (diffusion) time response were used in 
conjunction with the sums of squares approximation to get an estimate of the time 
constants which limit the overall response of both the PIN and the MQW structures. The 
diffusion time constant for the MQW structure (idiff- 1.5 ns) was found to be more than 
twice as large as that calculated for the PIN APD (T*ff~ 0.7 ns) and confirms that the 
undepleted MQW structure was diffusion-limited. 
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APPENDIX A 


AUGai. x As MATERIAL SYSTEM PARAMETERS 

GaAs/AlGaAs materials are frequently used in the fabrication of heterojunction 
devices. The parameters for GaAs can, in general, be derived from those for Al x Gai. x As 
by setting the value of x to zero. In the following few sections, the various optical and 
electrical parameters relevant to GaAs/AlGaAs structures modeling will be presented. 


A.1 A l QaAs Ba nd Pa r a m ete rs 

The net material bandgap is chosen as the minimum value of the various default bandgap 
valleys in AlGaAs given by the following equations : 

E g (G)=Ego + x(l . 1 55 + 0.37 x) [A- 1] 

E g (L)=1.734 + x(0.574 + 0.055 x) [A- 2] 

E g (X)=l .91 1 + x(0.005 + 0.245 x) [A- 3] 

The temperature dependence of the bandgap is introduced through the Ego parameter 
calculated from the following equation 1,2 : 


E g (T) = E g (0)~ 


aT 2 
T+p 


£,(300) + a 


300 2 
300 + /? 


T 2 

T + p 


[A- 4] 
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where, 

E,(0)= 1.519 eV 
a = 5.405 x 10-4 eV/K 
P = 204.0 K 

The AlGaAs electron affinity for a given x value can be calculated from: 

X = 4.07 - 0.85(£, (G) - E g ( GaAs )) [A- 5] 

and the density of states mass of the valence band is given by: 

m v = (0.082 3/2 + 0.4 5 3/2 ) 2/3 [A- 6] 

A. 2 AlGaAs Dielectric Properties 

The default value for the static dielectric constant for AlGaAs is given as a function of x 
by: 

s’ = 13.18 + 2.9 x [A- 7] 

In the case of GaAs, x=0, and e’ = £g,/J Eo = 13.18. 

A.3 AlGaAs Bandgap Narrowing 

Heavy doping of materials usually lead to bandgap narrowing which has an important 
effect on heterojunction characteristics. According to Lundstrom 3 , the bandgap narrowing 
effects are important only for p-type regions. The default bandgap narrowing parameters 
are shown in Table 1 for AlGaAs material. 
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Concentration 

cm"* 

Bandgap Narrowing 
m#V 

1.QE18 

31.0 

2.0E18 

36.0 

4.0E18 

44 2 

6.0E18 

48.6 

8.0E18 

61.7 

1.0E19 

64.3 

2.0E19 

61.1 

4.0E19 

64.4 

6.0E19 

61.9 

8.0E19 

66.9 

1.0E20 

63.2 

2.0E20 

18.0 


Table A-l : Default bandgap narrowing values 


A.4 AlGaAs Recombination Parameters 

The default parameters used for the various recombination mechanisms parameters are 
shown in Table A- 2 : 


Parameter 

Value 

Equation 

TnO 

l.OxlO' 9 

[2-19] 

TdO 

2.0X10- 8 

[2-201 


1.5xlO' 10 

[2-18] 

Cn 

5.0X10 30 

[2-17] 


l.OxlO" 31 

_[2 Z 17] 


Table A- 2: Default Recombination Parameters for AlGaAs 
A. 5 GaAs Impact Ionization Coefficients 

The coefficients used in the impact ionization model described by equations [2-21]-[2-23] 
are provided in Table A- 3. 
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Value 

■9 

1.899x10 s 

1 

5.75x10 s 

Pn 

1.82 

mm 

2.215x10 s 


6.57x10 s 

Pp 

1.75 


Table A- 3: Impact Ionization Coefficients for GaAs 
A. 5 AlGaAs Mobility Parameters 


The default low-field electron mobility for various ranges of Al x Gai.*As composition is 
given by the following set of equations : 


= 8000 - 1.818*10 4 x 
\i„ = 90 + 1 . 1 43 5 * 1 0 5 (x- 0.46) 2 
p„ = 90 + 3.75*10 4 (x- 0.46) 2 
H„ = 200-2. 0/(x-0. 46) 


(0< x < 0.429) 
(0.429 < x < 0.46) 
(0.46 <x< 0.5) 
(0.5 < x < 1.0) 


[A- 8] 
[A- 9] 
[A- 10] 
[A- 11] 


The default concentration-dependent mobility values for GaAs are given in Table A- 4. 


The following expression is used to implement a field-dependent mobility which provides a 


smooth transition between low-field and high-field behavior: 


= M> 


-|i <fi 


1 + 


MqE 


[A- 12] 
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Where 0 is a constant with a default value of 1.0 for both electrons and holes. The 



a 






































































APPENDIX B 


SAMPLE SIMULATION ALGORITHMS 
B .l Light IV Simulation of a PIN Structure 

S GaAs pin Simulation under Light (p=i=n=l micron) 

$by: Hicham Menkara 
S input deck : pin_bias6_lite.in 

$ 

S SECTION 1 : Mesh specification 

S 

mesh nx=3 ny=100 smooth=l diag.flip space.mult=1.0 master. out 

S 

Smesh reel nx=3 ny=400 smooth =4 diag.flip 


x.mesh i 

n=l loc= 

0 r=l 

x.mesh 

n=3 loc 

=75.0 

y.mesh 

n=l 

loc=0.0 

y.mesh 

n=8 

loc=0 96 

Sy.mesh 

n=45 

loc=1.2 

y.mesh 

n=50 

loc 3 15 

Sy.mesh 

n=255 

loc=3.3 

y.mesh 

n=92 

loc=2.04 

y.mesh 

n=100 

loc=3 


$ SECTION 2: Structure Specification 

$ 

region num=l GaAs y.max=l 
region num=2 GaAs y.min=l y.max=2 
region num=3 GaAs y.min=2 y.max=3 

S 

elec num=l x.min=0 x.max=7S y.max=0.0 
elec num=2 bot 
contact num=l con. resist 3 1 ©4 
contact num=2 conresist=le4 

$ 

doping region 3 1 uniform p.type conc=3e!8 
doping region 3 3 uniform n.type conc=3el8 

S 

$ SECTION 3: Material model specification 

$ 

material material=GaAs taupO=l e-9 taunO=l e-9 
impact material=GaAs selb 

models material=GaAs srh auger conmob fldmob print 

S 

$ SECTION 4: Optical source definition 

S 

beam num=l x.origin=37.5 y.origin=-1.0 angle 3 90.0 wavelength 3 . 6328 min.window=-9.0 max window=9.0 

S 
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$ SECTION S: Initial solution 

S 

symb newton carr=0 

solve init 

symb newton carr=2 

method trap autonr climit=75000 ctolt.fact=500.0 maxtrap=6 
solve prev 

Splot.2d x.min=0 x.max=75 y.min=0 y.max=3 grid depl.edge 
Splotld e.field a.x=37.5 b.x*37.5 a.y=0 b.y=5 points 
solve bl*0.001 

output e.field j.electron e. velocity e.mobility h.mobility qss c.temp h.temp val.band conband qfii qfp impact 
recomb tot.doping photogen flowlines u.auger u. radiative u.srh 
save outf=pin_nobias_lite7.out 

S 

S SECTION 6: Voltage ramp 

S 

log outf=pin_IV_lite7.1og master 

solve prev vl=0.0 vstep=-2 vfmal—10 elect=l 

solve project vl~12 vstep=-2 vfmal=-35 elect=l 

plot. Id alphan a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=eionl35.dat ascii 

plot. Id alphap a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=hionl35.dat ascii 

save outf==pin_bias35_lite7.out 

solve project vl— 35.2 vstep=-0.1 vfinal=-37.5 elect=l master 

solve project vl~37.6 vstep=-0.1 vfmal=-4 5 elect=l master 

plot. Id alphan a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=eionl38.dat ascii 

plot. Id alphap a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=hionl38.dat ascii 

save outf=pin_bias38_lite6.out 
tonyplot pin_IV_lite6.log 

S 

end 


B.2 CV Simulation of a doped-well MOW Structure 

S Doped MQW CV simulation example (10 barriers, 9 wells) 

S by Hicham Menkara 
S input deck : MQWCV7s2b.in 

S 

S SECTION 1: Mesh Specification 

S 

Smesh space. mult=4.0 

S 

mesh rect smooth =4 diag flip 
x mesh loc=0 s=75 

x. mesh loc=75.0 s=75 

y. mesh loc=0.0 s=0.4 

y.mesh loc=0.% s=0.4 
y.mesh loc=0.97 s=0.0024 
y.mesh loc=1.525 s=0.0024 

y.mesh loc=2.08 s=0.0024 
y.mesh loc=2.09 s=0.4 

y.mesh loc=3.05 s=0.4 

$ 

$ SECTION 2: Structure Specification 
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$ 

region number 5 1 x.min=0 x.max=75 y.min=0 y.max=l .0 gaas 
region number2 x.min=0 x.max=75 y.min= 1.0 y.max= 1.05 gaas 

region number 5 3 x.min=0 x.max=75 y.min=1.05 y.max=l.l material=AlGaAs x. composition =0.42 
region number=4 x.min=0 x.max=75 y.min=l.l y.max=1.15 gaas 

region number's x.min=0 x.max=75 y.min=1.15 y.max=1.2 material=AlGaAs x. com position =0.42 
region number=6 x.min=0 x.max=75 y.min=1.2 y.max=1.25 gaas 

region number 7 x.min=0 x.max=75 y.min=l .25 y.max= 1 .3 material=AJGaAs x.composition=0.42 
region number8 x.min=0 x.max=75 y.min=1.3 y.max=1.35 gaas 

region number9 x.min=0 x.max=75 y.min=1.35 y.max=1.4 material=AlGaAs x.composition=0.42 
region number 10 x.min=0 x.max=75 y.min=1.4 y.max=1.45 gaas 

region number 1 1 x.min=0 x.max=75 y.min=1.45 y.max=1.5 material=AlGaAs x. composition^. 42 
region number 12 x.min=0 x.max=75 y.min=1.5 y.max=1.55 gaas 

region number 13 x.min=0 x.tnax=75 y.min=1.55 y.max=1.6 material =AlGaAs xcomposition=0.42 
region number 14 x.min=0 x.max=75 y.min=1.6 y.max=1.65 gaas 

region number 15 x.min=0 x.max=75 y.min=1.65 y.max=1.7 material=AlGaAs x.composition=0.42 
region number 16 x.min=0 x.max=75 y.min=1.7 y.max=1.75 gaas 

region number 17 x.min=0 x.max=75 y.min=1.75 y.max=1.8 material=AlGaAs x.composition=0.42 
region number 18 x.min=0 x.max=75 y.min=1.8 y. mar 1.85 gaas 

region number 19 x.min=0 x.max=75 y.min=1.85 y. mar 1.9 material=AlGaAs x.composition=0.42 
region number 20 x.min=0 x.max=75 y min=1.9 y.max=195 gaas 

region number 21 x min=0 x.max=75 y.min=1.95 y.max=2.0 material=AlGaAs x. composition^. 42 
region number 22 x.min=0 x.max=75 y.min=2.0 y.max=2.05 gaas 
region number 2 3 x min=0 x.max=75 y.min=2.05 y.max=3.05 gaas 

* 

M #l=cathode #2=anode 
electrode name=cathode number 1 top 
electrode name= anode number 2 bottom 
« 

doping uniform conc=3e 18 p.type direction=y regions=l 
doping uniform conc=1.5e 18 p.type direction=y y.min=l. 115 y.max=l. 120 
doping uniform conc=1.5e 18 n.type direction=y y.min=1.130y.max=1.135 
doping uniform conc=1.5el8 p.type direction=y y.min=1.215 y.max=l. 220 
doping uniform conc=1.5el8 n.type direction=yy.min=l. 230 y.max=l. 235 
doping uniform conc=1.5el8 p.type direction=yy.min=l. 315 y.max=l. 320 
doping uniform conc=1.5el 8 n.type direction=y y.min= 1.330 y.max= 1.335 
doping uniform conc=1.5e 18 p.type direction=y y.min= 1.4 15 y.max= 1.420 
doping uniform conc=1.5e 18 n.type direction=yy.min= 1.430 y.max= 1.435 
doping uniform conc=1.5e 18 p.typedirection=yy.min=1.515 y.max=1.520 
doping uniform conc=1.5e 18 n.type direction=y y.min=l. 530 y.max=l. 535 
doping uniform conc=1.5e 18 p.type direction=y y.min=1.615 y.max=l. 620 
doping uniform conc=1.5el8 n.type directi on=y y.min= 1.630 y.max=1.635 
doping uniform conc=1.5el8 p.type direction=y y.min=1.715 y.max=1.720 
doping uniform conc=1.5el8 n.type direction=y y.min=1.730 y.max=1.735 
doping uniform conc=1.5e 18 p.type direction=y y.min=1.815 y.max=l. 820 
doping uniform conc=1.5e 18 n.type direction=y y.min=l. 830 y.max=l. 835 
doping uniform conc=1.5el8 p.type direction=y y. min= 1.91 5 y.max= 1.920 
doping uniform conc=1.5el8 n.type direction=y y.min=1.930 y.max=l. 935 
doping uniform conc=3e 18 n.type direction=y regions=23 

$ 

$ SECTION 3: Material Model Specification 

S 

material taup0=2.e-6 taun0=2.e-6 
models srh auger conmob fldmob 
impact selb 
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I 

S SECTION 4: Initial Solution 

S 

symb Newton carr=2 
method comb trap autonr 
solve init 

output e.field j. electron e. velocity e.mobility h.mobility qss e.temp h.temp val.band con. band qfh qfp impact 

recomb tot doping 

save outf=nobiascent7s2b.out 

log outfile=MQWC V3j2b.log master 

method itlimit*50 autonr nrcriterion=0.1 trap atrap“0.5 maxtrap=10 

solve vl=0 v2=0 vstep=-0 1 vfinal=-10 elect=l ac direct freq=le6 vss*0.030 terminal=l 
solve vl*-ll v2=0 vstep=-l vftnal=-25 elect* 1 ac direct freq=lc6 vss=0.030 terminal=l 
save outf=MQWcent725s2b.out 
tonyplot MQWCV3j2b.log 
end 


B.3 Spectral Response Simulation of a PIN Structure 

S Spectral response of PIN structure (p=i=n= 1 micron) 

S by Hicham Menkara 
S input deck : PINspec2a.in 

S 

$ SECTION 1: Mesh Specification 

S 

Smesh space. mult=4.0 

S 

mesh rect smooth =4 diag.flip 
x.mesh loc=0 s=75 

x. mesh loc=75.0 s=75 

y. mesh loc=0.0 s=0.01 

y.mesh loc=0.96 s=0.01 
y.mesh loc=0.97 s=0.004 
y.mesh loc=1.05 s=0.01 
y.mesh loc=2.08 s=0.004 
y.mesh loc=2.09 s=0.4 

y.mesh loc=3.05 s=0.4 
$ 

S SECTION 2: Structure Specification 

$ 

region number* 1 x.min=0 x.max=75 y.min*0 y.max* 1 .0 gaas 
region number=2 x.min=0 x.max=75 y.min= 1 .0 y.max=2.05 gaas 
region number=3 x.min=0 x.max=75 y.min=2.05 y.max=3.05 gaas 

# 

# #l=cathode #2= anode 

electrode name*cathode number* 1 top 

electrode name=anode number=2 bottom 

doping uniform conc=3el8 p.type direction*y regions* 1 
doping uniform conc=3el8 n. type direction=y regions* 3 

S 

S SECTION 3: Material Model Specification 

S 
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material material=AlGaAs taupO=2.e-8 taunO=l.e-9 copt=1.5e-10 augn=5e-30 augp=le-31 

impact aelber anl=1.899e5 an2=1.899e5 bnl=5.75e5 bn2=5.75e5 apl=221500 ap2=221500 bpl=657000 

bp2=657000 betan=1.82 betap=1.75 egran=0 

models material=GaAs srh auger conmob fldmob print 

S 

S SECTION 4: Optical source definition 

$ 

beam num=l x.origin=37.5 y. origin =-1.0 angle=90.0 wavelength=.2 min.window=-9.0 max window=9.0 

S 

$ SECTION 5: Initial solution 

S 

symb newton can=0 

solve init 

symb newton can=2 

method trap autonr climit=75000 ctolt.fact=500.0 maxtrap= 1 0 
solve prev 

output e. field j.electron e. velocity e.mobility h.mobility qss e.temp h.temp val.band con. band qfn qfp impact 
recomb tot.doping photogen traps flowlines ey. velocity ex. velocity hx. velocity hy. velocity u.auger u.srh u.radiative 
solve prevbl=l 

save outf=PINspec2a0.out 

S 

S SECTION 6: spectral response 

S 

log outf=PINspec2a.log 

solve prev bl=l lambda=0.2 

solve prev bl=l lambda=0,225 

solve prev bl— 1 lambda=0.25 

solve prev bl=l lambda=0.275 

solve prev bl = 1 lambda=0.3 

solve prev bl=l lambda=0.325 

solve prev bl=l lambda=0,35 

solve prev bl = 1 lambda=0.375 

solve prev bl=l lambda=0.4 

save outf=PINspec400.out 
solve prev bl=l lambda=0.425 

solve prev bl = 1 lambda =0 4 5 

solve prev bl = 1 lambda=0 475 

solve prev bl =1 lambda=0.5 

solve prev bl = 1 lambda=0.525 

solve prev bl = 1 lambda=0.55 

solve prev bl=l lambda=0.575 

solve prev bl =1 lambda=0.6 

solve prev b 1 = 1 lambda=0.625 

solve prev bl =1 lambda=0.6325 

save outf=PINspec632.out 
solve prev bl =1 lambda=0.65 

solve prev bl = 1 lambda=0.675 

solve prev bl=l lambda=0.7 

solve prev bl=l lambda=0.725 

solve prev bl =1 lambda=0.75 

solve prev bl=l lambda =0.77 5 

solve prev bl=l lambda=0.8 

solve prev bl=l lambda=0 825 

solve prev b 1 = 1 lambda=0.85 

save outf=PINspec850.out 
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solve prev bl=l lambda=0.865 

solve prev bl=l lambda =0 875 

solve prev bl=l Lambda=0 885 

solve prev bl=l lambda =0 9 

solve prev bl=l lambda =0 91 

solve prev bl=l lambda=0.9184 

save outf=PINspec2a2.out 
tonyplot PINspec2a.log 
end 
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APPENDIX A 


Al x Gai.„ As MATERIAL SYSTEM PARAMETERS 

GaAs/AlGaAs materials are frequently used in the fabrication of heterojunction 
devices. The parameters for GaAs can, in general, be derived from those for Al x Gai_ x As 
by setting the value of x to zero. In the following few sections, the various optical and 
electrical parameters relevant to GaAs/AlGaAs structures modeling are presented. 


A. 1 AlGaAs Band Parameters 

The net material bandgap is chosen as the minimum value of the various default 
bandgap valleys in AlGaAs given by the following equations : 

E g (G)=Ego + x(l . 1 55 + 0.37 x) [A- 1] 

E,(L)=1 .734 + x(0.574 + 0.055 x) [A- 2] 

E g (X)=1.91 1 + x(0.005 + 0.245 x) [A- 3] 

The temperature dependence of the bandgap is introduced through the Ego parameter 
calculated from the following equation 56,57 : 


aT 2 


E g (T) = E t { 0) - — = £,(300) + a 


300 2 T 2 
300 + /? ~ T + 0 


[A- 4] 
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where, 

E,(0)= 1.519 eV 
a = 5.405 x 10-4 eV/K 
P = 204.0 K 

The AlGaAs electron affinity for a given x value can be calculated from: 

X = 4.07 - 0.85(£, (G) - E g (GaAs)) [A- 5] 

and the density of states mass of the valence band is given by: 

m v = (0.082 3/2 + 0.45 3 ' 2 ) 2/3 [A- 6] 

A. 2 AlGaAs Dielectric Properties 

The default value for the static dielectric constant for AlGaAs is given as a 
function of x by: 

s' - 13.18 + 2.9 x [A- 7] 

In the case of GaAs, x=0, and e' = Zg,\J £o = 13.18. 

A.3 AlGaAs Bandgap Narrowing 

Heavy doping of materials usually lead to bandgap narrowing which has an 
important effect on heterojunction characteristics. According to Lundstrom 51 , the bandgap 
narrowing effects are important only for p-type regions. The default bandgap narrowing 
parameters are shown in Table 1 for AlGaAs material: 
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Table A-l : Default bandgap narrowing values 


Concentration 

cm 

Bandgap Narrowing 
maV 

L0E18 

31.0 

2.0E18 

38.0 

4.0E18 

44.2 

6.0E18 

48.5 

8.0E18 

51.7 

L0E19 

54.3 

2.0E19 

61.1 

4.0E19 


6.0E19 

61.9 

8.0E19 

66.9 

1.0E20 

53.2 

2.0E20 

18.0 


A.4 AlGaAs Recombination Parameters 

The default parameters used for the various recombination mechanisms parameters are 
shown in Table A-2 : 


Table A-2: Default Recombination Parameters for AlGaAs 


Parameter 

Value 

Equation 

tnO 

l.OxlO" 9 

12-19] 

XpO 

2.0x10-® 

12-201 

Co* 

1.5xl0 10 

12-181 

Cn 

5.0X10 30 

12-17] 

Cp 

l.OxlO" 3 ' 



A. 5 GaAs Impact Ionization Coefficients 

The coefficients used in the impact ionization model described by equations [2- 
21]-[2-23] are provided in Table A-3. 

Table A-3: Impact Ionization Coefficients for GaAs 


134 



























mm 

Value 

mm 

1.899xl0 5 


5.75xl0 5 

A 

1.82 

mm 

2.215x10 s 

■n 

6.57x10 s 

A 

1.75 


A. 6 AlGaAs Mobility Parameters 

The default low-field electron mobility for various ranges of Al x Gai. x As composition is 


given by the following set of equations : 


p„ = 8000 - 1.818*10 4 x 
p„ = 90 + 1.1435*10 5 (x- 0.46) 2 
|i„ = 90 + 3.75*10 4 (x-0.46) 2 
V* = 200-2.0/(x-0.46) 


(0< x < 0.429) [A- 8] 

(0.429 < x < 0.46) [A- 9] 

(0.46 < x < 0.5) [A- 10) 

(0.5 < x < 1.0) [A- 11] 


The default concentration-dependent mobility values for GaAs are given in Table A-4. 


The following expression is used to implement a field-dependent mobility which provides a 
smooth transition between low-field and high-field behavior: 


M(E) = Mi 




1 + 






[A- 12] 
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Where P is a constant with a default value of 1 .0 for both electrons and holes. The 


saturation velocities for AlGaAs are given by: 

= 1.13* 10 7 - 1.2* 10 4 T [A- 13] 

= 1.12* 10 7 - 1.2* 10 4 T [A- 14] 


Table A-4: Default concentration-dependent mobilities for GaAs 



Mobility In GaAs 
(cm*/v-8) 

Conctntntion (cm-3) 

Electrons 

HoU» 

1 . 0 * 10 14 

8000.0 

390.0 

2 . 0 * 1 0 1 4 

7718.0 

380.0 

4 . 0 * 10' 4 

7445.0 

375.0 

6 . 0 * 10 14 

7290.0 

360.0 

8 . 0 * 10 14 

7182.0 

350.0 

1 . 0 * I 0 1S 

7300.0 

340.0 

2 . 0 * 10 15 

6847.0 

335.0 

4 . 0*1 O ' 5 

6422.0 

320.0 

6 . 0 * 1 0 1 5 

6185.0 

315.0 

80 * 10 15 

6023.0 

305.0 

1 . 0 * I 0 IS 

5900.0 

302.0 

2 . 0 * I O 14 

5474.0 

300.0 

4 . 0*1 O ' 6 

5079.0 

285.0 

6 . 0 * 10 " 

4861.0 

270.0 

8 . 0 * 1 0 1 * 

4712.0 

245.0 

1 . 0 * 10 n 

4600.0 

240.0 

2 . 0 * 10 17 

3874.0 

210.0 

4 . 0 * 10 17 

3263.0 

205 0 

6 . 0 * 10 17 

2950.0 

200 0 

8 . 0 * 10' 7 

2747.0 

1869 

1 . 0 * 10 " 

2600.0 

170.0 

2 . 0 * 10 " 

2060.0 

130 0 

4 . 0 * 10 " 

1632.0 

90.0 

6 . 0 * 10 " 

14240 

74.5 
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APPENDIX B 


SAMPLE SIMULATION ALGORITHMS 

B .l Light IV Simulation of a PIN Structure 

S GaAs pin Simulation under Light (p= c i*n=l micron) 

Sby: Hicham Menkara 
S input deck : pin_bias6_lite.in 

S 

S SECTION 1 : Mesh specification 

S 

mesh nx=3 ny=100 smooth=l diag.flip space. mult= 1.0 master.out 

S 

Smesh rect nx=3 ny=400 smooth=4 diag.flip 
x.mesh n=l loc=0 r=l 


x.mesh 

n=3 loc=75.0 

y.mesh 

n=l 

loc=0.0 

y.mesh 

n=8 

loc=0.96 

Sy.mesh 

n*45 

loc=1.2 

y.mesh 

n=50 

loc=1.5 

Sy.mesh 

n=255 

loc=3.3 

y.mesh 

n=92 

loc=2.04 

y.mesh 

n=100 

loc=3 


S SECTION 2: Structure Specification 

$ 

region num=l GaAs y.max=l 
region num=2 GaAs y.min=l y.max=2 
region num=3 GaAs y.min=2 y.max=3 

S 

elec num=l x.min=0 x.max=75 y.max=0.0 
elec num=2 bot 
contact num=l con.resist=le4 
contact num=2 con.resist=le4 

S 

doping regional uniform p.type conc=3el 8 
doping region^ 3 uniform n.type conc=3el8 

S 

$ SECTION 3: Material model specification 

S 

material material=GaAs taupfr* 1 .e-9 taunO* 1 e-9 
impact material=GaAs selb 

models material=GaAs srh auger conmob fldmob print 

$ 

$ SECTION 4: Optical source definition 

S 
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beam num=l x.origin=37.5 y. origin— 1.0 angle=90.0 wavelength 2 . 6328 min . window=-9 . 0 max.window=9.0 

$ 

$ SECTION 5: Initial solution 

S 

symb newton carrO 

solve init 

symb newton carr=2 

method trap autonr climit=75000 ctolt.fact=500.0 maxtrap=6 
solve prev 

Splot2d x.min=0 x.max=75 y,min=0 y.max=3 grid depl.edge 
Splot.ld e. field a.x=37.5 b.x=37.5 a.y=0 b.y=5 points 
solve bl =0.001 

output e.field j. electron e. velocity e.mobility h. mobility qss e.temp h.temp val.band con. band qfh qfp impact 
recomb tot. doping photogen flowlines u.auger u. radiative u.srh 
save outf=pin nobias_lite7 out 

S 

$ SECTION 6: Voltage ramp 

S 

log outf=pin_IV_lite7.1og master 

solve prev vl=0.0 vstep=-2 vfinal— 10 elect 2 1 

solve project v 1—1 2 vstep=-2 vfinal=-35 elect 2 1 

plot. Id alphan a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=eion!35.dat ascii 

plot. Id alphap a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=hionl35.dat ascii 

save outf==pin_bias35_lite7.out 

solve project vl— 35.2 vstep=-0.1 vfinal— 37.5 elect 2 1 master 

solve project vl— 37,6 vstep=-0.1 vfinal=-4 5 elect 2 1 master 

plot. Id alphan a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=eionl38.dat ascii 

plot ld alphap a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=hionl38.dat ascii 

save outf=pin_bias38_lite6.out 
tonyplot pin_IV_lite6.log 

S 

end 


B.2 CV Simulation of a doped-well MOW Structure 

S Doped MQW CV simulation example (10 barriers, 9 wells) 

S by Hicham Menkara 
S input deck : MQWCV7s2b.in 

S 

$ SECTION 1: Mesh Specification 

S 

Smesh space. mult=4.0 

S 

mesh rect smooth =4 diag.flip 
x.mesh loc=0 s=75 

x. mesh loc=75.0 s=75 

y. mesh loc=0.0 s=0.4 

y.mesh loc=0.96 s=0.4 
y.mesh loc=0.97 s=0.0024 
y.mesh loc=1.525 s=0.0024 

y.mesh loc=2 08 s=0.0024 
y.mesh loc=2.09 s=0.4 

y.mesh loc=3.05 s=0.4 
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s 

S SECTION 2: Structure Specification 

$ 

region number=l x.min=0 x.max=75 y.min=0 y.max=1.0 gaas 
region number=2 x.min=0 x.max=75 y.min= 1 .0 y.max=l .05 gaas 

region number=3 x.min=0 x.max=75 y.min* 1 .05 y.max=l. 1 mated al=AJGaAs x.composition=0.42 
region number=4 x.min=0 x.max=75 y.min* 1.1 y.max=1.15 gaas 

region number* 5 x.min=0 x.max=75 y.min* 1.1 5 y.max=1.2 material=AlGaAs x.composition=0.42 
region number=6 x.min=0 x.max*75 y.min* 1.2 y.max* 1.25 gaas 

region number*7 x.min=0 x.max*75 y.min=1.25 y. max* 1.3 material*AlGaAs x.composition=0 42 
region number*8 x min=0 x.tnax=75 y.min* 1.3 y.max* 1.3 5 gaas 

region number*9 x.min=0 x.max=75 y.min* 1 .35 y.max*1.4 material*AlGaAs x.composition=0.42 
region number* 10 x.min=0 x.max=75 y.min* 1.4 y. max* 1.45 gaas 

region number* 1 1 x.min*0 x.max*75 y.min=l .45 y. max* 1.5 material*AlGaAs xcomposition=0.42 
region number* 12 x.min=0 x.max=75 y.min* 1.5 y.max=1.55 gaas 

region number* 13 x.min=0 x.max=75 y.min=1.55 y. max* 1.6 material=AlGaAs x.composition=0.42 
region number* 14 x.min=0 x.max=75 y.min* 1.6 y. max* 1.65 gaas 

region number* 15 x.min=0 x.max=75 y.min* 1.65 y. max* 1.7 material=AlGaAs x. composition^. 42 
region number* 16 x.min=0 x.max=75 y.min* 1.7 y.max=1.75 gaas 

region number* 17 x.min=0 x.max=75 y.min* 1.75 y. max* 1.8 material* AlGaAs x. composition^. 42 
region number* 18 x.min=0 x.max=75 y.min* 1.8 y.max=1.85 gaas 

region number* 19 x.min=0 x.max*75 y.min* 1.85 y. max* 1.9 material*AlGaAs x.composition=0.42 
region number=20 x.min=0 x.max*75 y.min* 1.9 y.max=1.95 gaas 

region number=21 x.min=0 x.max=75 y.min* 1.95 y.max=2 0 material=AlGaAs x.composition=0.42 
region number=22 x.min=0 x.max=75 y.min=2.0 y.max=2.05 gaas 
region number=23 x.min=0 x.max=75 y.min=2.05 y.max=3.05 gaas 

« 

# #l=cathode #2*anode 

electrode name* cathode number* 1 top 
electrode name=anode number* 2 bottom 

# 

doping uniform conc=3e 18 p .type directi on =y regions* 1 
doping uniform cone* 1 . 5e 1 8 p.type direction=y y.min* 1.1 15 y.max=1.120 
doping uniform cone* 1.5e 18 n.typedirection=y y.min=1.130y.max=1.135 
doping uniform cone* 1 . 5e 1 8 p.type direction=y y.min=1.215 y.max=1.220 
doping uniform cone* 1.5e 18 n.type direction=y y.min* 1.230 y.max=1.235 
doping uniform conc=1.5el8 p.type direction=y y.min=l. 315 y.max=l. 320 
doping uniform cone* 1.5e 18 n.type direction=y y.min* 1. 330 y.max=l. 335 
doping uniform conc=1.5el8 p.type direction=y y.min* 1.4 15 y. max* 1.420 
doping uniform cone* 1.5e 18 n.type direction=y y.min* 1.430 y. max* 1.435 
doping uniform cone* 1.5e 18 p.type direction=yy.min=1.515y.max=l. 520 
doping uniform conc=1.5el8 n.type direction=y y.min* 1.5 30 y. max* 1.535 
doping uniform cone* 1 . 5e 1 8 p.type direction=y y.min=1.615 y.max=1.620 
doping uniform cone* 1.5e 18 n.type direction=y y.min* 1 .630 y.max=1.635 
doping uniform cone* 1.5c 18 p.type direction=y y.min=1.715 y.max=1.720 
doping uniform conc=1.5el8 n.type direction=yy.min=1.730y.max=l. 735 
doping uniform cone* 1.5c 18 p.type direction*y y.min=1.815 y. max* 1.820 
doping uniform conc=15el8 n.type direction=y y.min* 1. 830 y.max=l. 835 
doping uniform conc=1.5el8 p.type dircction=y y.min* 1.9 15 y. max* 1.920 
doping uniform cone* 1.5e 18 n.type direction=y y.min* 1.930 y. max* 1.935 
doping uniform conc*3el 8 n.type direction=y regions=23 

S 

$ SECTION 3: Material Model Specification 

$ 

material taup0=2.e-6 taun0=2.e-6 
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models srh auger conmob fldmob 
impact selb 
$ 

$ SECTION 4: Initial Solution 

S 

symb Newton cair=2 
method comb trap autonr 
solve init 

output e.field j. electron e. velocity e.mobility h.mobility qss e.temp h.temp val.band con. band qfe qfp impact 
recomb toLdoping 

save outf=nobiascent7s2b.out 

log outfile=MQWCV3j2b.log master 

method itlimit=50 autonr nrcriterion=0.1 trap atrap=0.5 maxtrap=10 

solve vl=0 v2=0 vstep=-0.1 vfinal=-10 elect=l ac direct freq=le6 vss=0 .030 terminal=l 
solve vl=-ll v2=0 vstep=-l vfinal=-25 elect=l ac direct freq= 1 e6 vss=0.030 terminal=l 
save outf=MQWcent725s2b.out 
tonyplot MQWCV3j2b.log 
end 


B,3 Spectral Response Simulation of a PIN Structure 

SSpectral response of PIN structure (p=i=n= 1 micron) 

S by Hicham Menkara 
Sinput deck : PINspec2a.in 

S 

S SECTION 1: Mesh Specification 

S 

Smesh space.mult=4.0 

S 

mesh rect smooth =4 diag.flip 
x.mesh loc=0 s=75 

x. mesh loc=75.0 s=75 

y. mesh locO.O s=0.01 

y.mesh loc=0.96 s=0.01 
y.mesh loc=0 97 s=0.004 
y.mesh loc=1.05 s=0.01 
y.mesh loc=2.08 s=0.004 
y.mesh loc=2.09 s=0.4 

y.mesh loc=3.05 s=0.4 
S 

$ SECTION 2: Structure Specification 

S 

region number=l x.min=0 x.max=75 y.min=0 y.max=1.0 gaas 
region number=2 x.min=0 x.max=75 y.min=l .0 y.max=2.05 gaas 
region number=3 x.min=0 x.max=75 y.min=2.05 y.max=3.05 gaas 

* 

# #l=cathode #2=anode 

electrode name=cathode number=l top 

electrode name=anode number=2 bottom 

doping uniform conc=3e 18 p.type direction=y regions=l 
doping uniform conc=3el8 n.type direction=y regions=3 

S 
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S SECTION 3: Material Model Specification 

S 

material material AlGaAs taupO=2.e-8taunO=l.e-9 copt=1.5e-10 augn=5e-30 augp=le-31 

impact aelber anl=1.899e5 sn2=1.899e5 bnl=5.75e5 bn2=5.75e5 apl=221500 ap2=221500 bp 1=657000 

bp2=657000 betan=1.82 betap=1.75 egran=0 

models materia] =GaAs srh auger conmob fldmob print 

S 

S SECTION 4: Optical source definition 

S 

beam num=l x.origin=37.5 y.origin=-1.0 angle=90.0 wavelength=.2 min window=-9 0 max window=9.0 

S 

S SECTION 5: Initial solution 

S 

symb newton carr=0 

solve init 

symb newton carr=2 

method trap autonr climit=75000 ctolt.fact=500.0 maxtrap=10 
solve prev 

output e.field j. electron e. velocity e.mobility h.mobility qss e.temp h.temp val.band con. band qfii qfp impact 
recomb tot.doping photogen traps flowlines ey. velocity ex. velocity hx. velocity hy. velocity u.auger u.srh u.radiative 
solve prevbl=l 

save outf=PINspec2a0.out 

S 

$ SECTION 6: spectral response 

S 

log outfHPINspec2a.log 

solve prev bl=l lambda=0.2 

solve prev bl=l Lambda =0 22 5 

solve prev bl=l lambda =0.25 

solve prev bl=l lambda =0.2 75 

solve prev bl=l lambda=0.3 

solve prev bl=l lambda=0.325 

solve prev bl=l lambda=0.35 

solve prev bl=l lambda=0.375 

solve prev bl=l lambda=0 4 

save outf=PINspec400.out 
solve prev bl = 1 lambda=0.425 

solve prev bl =1 lambda=0 45 

solve prev bl = 1 lambda=0.475 

solve prev bl =1 lambda =0.5 

solve prev bl=l lambda=0.525 

solve prev bl=l lambda=0.55 

solve prev bl = 1 lambda=0.575 

solve prev bl=l lambda=0.6 

solve prev bl=l lambda=0.625 

solve prev bl=l lam bda=0. 6325 

save outf=PINspec632.out 
solve prev bl = 1 lambda=0.65 

solve prev bl=l lambda=0.675 

solve prev bl = 1 lambda=0.7 

solve prev bl = 1 lambda=0.725 

solve prev bl = 1 lambda=0.75 

solve prev bl = 1 lambda=0.775 

solve prev bl=l lambda=0.8 

solve prev bl=l lambda=0.825 
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solve prev bl=l lambda =0.8 5 
save outf=PINspec8 50. out 
solve prev bl = 1 lambda =0 86 5 

solve prev bl = 1 lambda=0 875 

solve prev bl=l lambda=0.885 

solve prev bl=l lambda =0 9 

solve prev bl=l lambda =0.91 

solve prev bl=l lambda=0.9184 

save outf=PIN spec2a2 . out 
tonyplot PINspec2a.log 
end 
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APPENDIX B 


SAMPLE SIMULATION ALGORITHMS 
B. 1 Light IV Simulation of a PIN Structure 

$ GaAs pin Simulation under Light (p=i=n=l micron) 

Sby: Hicham Menkara 
$ input deck : pin_bias6_lite.in 

S 

S SECTION 1: Mesh specification 

$ 

meshnx=3 ny=100 smooth=l diag.flip space. mult= 1.0 master out 

S 

Smesh rect nx=3 ny=400 smooth =4 diag.flip 
x.mesh n=l loc=0 r=l 


x.mesh 

n=3 loc= 

=75.0 

y.mesh 

n=l 

loc=0.0 

y.mesh 

n=8 

loc=0.96 

Sy.mesh 

n=45 

loc=1.2 

y.mesh 

n=50 

loc=1.5 

Sy.mesh 

n=255 

loc=3.3 

y.mesh 

n=92 

loc=2.04 

y.mesh 

n=100 

loc=3 


S SECTION 2: Structure Specification 

S 

region num=l GaAs y.max=l 
region num=2 GaAs y.min=l y.max=2 
region num=3 GaAs y.min=2 y.max=3 

S 

elec num=l x.min=0 x.max=75 y.max=0.0 
elec num=2 bot 
contact num=l con.resist=le4 
contact num=2 con.resist=le4 

S 

doping region=l uniform p.type conc=3el8 
doping region=3 uniform n.type conc=3el 8 

S 

S SECTION 3: Material model specification 

S 

material material=GaAs taup(M.e-9 taun0=l.e-9 
impact material-GaAs selb 

models material=GaAs srh auger conmob fldmob print 

S 

$ SECTION 4: Optical source definition 

$ 

beam num=l x.origin=37.5 y. origin—]. 0 angle=90.0 wavelength=.6328 min.window=-9.0 max.window=9.0 

$ 



S SECTION 5: Initial solution 

S 

symb newton carr=0 

solve init 

symb newton carr=2 

method trap autonr climit=7S000 ctolt.fact= 500.0 maxtrap=6 
solve prev 

Splot2d x.min=0 x.max=75 y.min=0 y.max=3 grid depl.edge 
Splotld e.field a.x=37.5 b.x=37.5 a.y=0 b.y=5 points 
solve bl =0.001 

output e.field j.electron e. velocity e.mobility h.mobility qss e.temp h.temp val.band con. band qfii qfp impact 
recomb tot doping photogen flowlines u. auger u. radiative u.srh 
save outf=pin nobias_lite7.out 

S 

S SECTION 6: Voltage ramp 

$ 

log outf=p in_IV_lite7.log master 

solve prev vl=0.0 vstep=-2 vfinal=-10 elect=l 

solve project vl =-12 vstep=-2 vfinal=-35 elect=l 

plot. Id alphan a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=eionl35.dat ascii 

plot.ld alphap a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=hionl35.dat ascii 

save outf==pin_bias35_lite7,out 

solve project vl=-35.2 vstep=-0.1 vfinal=-37. 5 elect=l master 

solve project vl— 37.6 vstep=-0.1 vfinal=-45 elect=l master 

plot. Id alphan a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile=eionl38.dat ascii 

plotld alphap a.x=37.5 b.x=37.5 a.y=0 b.y=5 points outfile= hionl38.dat ascii 

save outf=pin_bias38_lite6.out 
tonyplot pin_IV_lite6.log 

$ 

end 


B.2 CV Simulation of a doped-well MOW Structure 

S Doped MQW CV simulation example (10 barriers, 9 wells) 

$ by Hicham Menkara 
S input deck : MQWCV7s2b.in 

S 

S SECTION 1: Mesh Specification 

S 

Smesh space. mult=4.0 

S 

mesh rect smooth =4 diag flip 
x.mesh loc=0 s=75 

x. mesh loc=75.0 s=75 

y. mesh loc=0.0 s=0.4 

y.mesh loc=0.96 s=0.4 
y.mesh loc=0.97 s=0.0024 
y.mesh loc= 1.525 s=0.0024 

y.mesh loc=2.08 s=0.0024 
y.mesh loc=2 09 s=0.4 

y.mesh loc=3.05 s=0.4 



$ 

$ SECTION 2: Structure Specification 

S 

region number* 1 x.min=0 x.max=75 y.min*0 y. max* 1.0 gaas 
region number*2 x.min=0 x.max*75 y.min* 1.0 y.max* 1.05 gaas 

region number* 3 x.min=0 x.max=75 y.min=l .05 y.max* 1 . 1 material*AlGaAs x.composition=0.42 
region number=4 x.min=0 x.max=75 y.min* 1.1 y.max* 1.1 5 gaas 

region number* 5 x.min*0 x.max=75 y.min* 1.15 y.max= 1 .2 material = A1 GaAs x,composition=0 .42 
region number*6 x.min*0 x.max=75 y.min=1.2 y. max* 1.25 gaas 

region number=7 x.min=0 x.max=75 y.min=l .25 y.max=l .3 material=AlGaAs x.composition=0.42 
region number=8 x.min*0 x.max=75 y.min=1.3 y.max=1.35 gaas 

region number* 9 x.min=0 x.max=75 y.min=l .35 y.max=l .4 material = AlGaAs x.composition=0.42 
region number* 10 x.min=0 x.max=75 y.min* 1.4 y. max* 1.4 5 gaas 

region number* 1 1 x.min=0 x.max=75 y. min* 1.45 y. max* 1.5 material* A1 GaAs x. composition^. 4 2 
region number* 12 x.min=0 x.max=75 y.min* 1.5 y. max* 1.55 gaas 

region number* 13 x min=0 x.max*75 y.min* 1.55 y.max* 1.6 material* AlGaAs x.composition=0.42 
region number* 14 x.min=0 x.max=75 y.min* 1.6 y. max* 1.65 gaas 

region number* 1 5 x min=0 x.max=75 y.min* 1 .65 y .max* 1 .7 material* AlGaAs x.composition=0.42 
region number* 16 x.min=0 x.max=75 y.min* 1.7 y.max=1.75 gaas 

region number* 17 x.min=0 x.max=75 y.min* 1.75 y. max* 1.8 material=AlGaAs x composition*0.42 
region number* 18 x.min=0 x.max=75 y.min* 1.8 y. max* 1.85 gaas 

region number* 19 x min*0 x max=75 y.min* 1.85 y.max=1.9 material* AlGaAs x. composition =0.42 
region number=20 x.min=0 x.max=75 y. min* 1. 9 y. max* 1.95 gaas 

region number=21 x.min=0 x.max*75 y.min* 1.95 y.max=2.0 material=A)GaAs x.composition=0.42 
region number=22 x.min=0 x.max=75 y.min=2.0 y.max=2.05 gaas 
region number=23 x.min=0 x.max=75 y.min=2.05 y.max=3.05 gaas 
M 

# #l=cathode #2=anode 

electrode name=cathode number* 1 top 
electrode name* an ode number* 2 bottom 

# 


doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 

doping 


uniform cone* 3e 18 p 
uniform conc=1.5el8 
uniform conc=1.5el 8 
uniform conc=1.5el8 
uniform conc*l .5el 8 
uniform conc=1.5el 8 
uniform conc=1.5el8 
uniform conc=1.5el8 
uniform conc=1.5el8 
uniform conc*1.5el8 
uniform «mc=1.5el8 
uniform conc=1.5el8 
uniform conc=1.5el8 
uniform conc*1.5el8 
uniform conc=1.5el8 
uniform conc*1.5el8 
uniform conc*1.5el8 
uniform conc*1.5el8 
uniform conc=1.5el8 
uniform conc=3el8 n. 


type direction =y regions* 1 
p.type direction=y y.min* 1.115 
n.tyjx: direction^ y.min* 1.130 
p.type direction=y y.min=1.215 
n.type direction=y y.min* 1.230 
p.type direction =y y.min=1.315 
n.type directi on=y y.min* 1.330 
p.type direction=y y.min* 1.4 15 
n.type direction=y y.min* 1.430 
p.type direction=y y.min=l .5 1 5 
n.type direction*y y.min* 1.530 
p.type direction*y y.min* 1.6 15 
n.type direction=y y.min* 1.630 
p.type direction=y y.min=1.715 
n.type directi on =y y.min* 1.730 
p.type direction=y y.min* 1.8 15 
n.type directi on=y y.min=1.830 
p.type direction=y y.min=1.915 
n.type direction*y y.min=1.930 
type direction *y regions=23 


y.max* 1.120 
y.max=1.135 
y.max* 1.220 
y.max=1.235 
y.max=1.320 
y.max=1.335 
y.max* 1.420 
y.max=1.435 
y.max* 1.520 
y.max* 1.535 
y.max* 1.620 
y.max* 1.635 
y.max* 1.720 
y.max* 1.735 
y.max=1.820 
y.max*1.835 
y.max=1.920 
y.max=1.935 


{ 

S SECTION 3: Material Model Specification 

S 

material taup0=2.e-<> taun0=2.e-6 
models srh auger conmob fldmob 
impact selb 
S 

S SECTION 4: Initial Solution 
S 

symb Newton carr=2 
method comb trap autonr 
solve init 

output e.field j. electron e. velocity e.mobility h.mobility qss e.temp h.temp val.band conband qfh qfp impact 
recomb tot doping 

save outf=nobiascent7s2b.out 

log outfile=MQWC V3j2b.log master 

method itlimit=50 autonr nrcriterion=0.1 trap atrap=0.5 maxtrap=10 

solve vl=0 v2=0 vstep=-0.1 vfinal=-10 elect=l ac direct freq=le6 vss=0.030 terminal 1 
solve vl— 11 v2=0 vstep=-l vfmal=-25 elect=l ac direct freq=le6 vss=0.030 terminal 1 
save outf=MQWcent725s2b.out 
tonyplot MQWCV3j2b.log 
end 


B.3 Spectral Response Simulation of a PIN Structure 

{Spectral response of PIN structure (p=i=n= 1 micron) 

S by Hicham Menkara 
{input deck : PINspec2a.in 

$ 

$ SECTION 1: Mesh Specification 

S 

{mesh space.mult=4.0 

{ 

mesh rect smooth =4 diag flip 
x.mesh loc=0 s=75 

x. mesh loc=75.0 s=75 

y. mesh loc=0.0 s=0 01 

y.mesh loc=0.96 s=0 01 
y.mesh loc=0.97 s=0.004 
y.mesh loc=1.05 s=0.01 
y.mesh loc=2.08 s=0.004 
y.mesh loc=2.09 s=0.4 

y.mesh loc=3.05 s=0.4 

{ 

{ SECTION 2: Structure Specification 

{ 

region number^ 1 x.min=0 x.max=75 y.min=0 y.max= 1 .0 gaas 
region number=2 x.min=0 x.max=75 y,tnin=1.0 y.max=2.05 gaas 
region number=3 x.min=0 x.max=75 y.min=2.05 y.max=3.05 gaas 

# 



# #l=cathode #2=anode 

electrode name=cathode number=l top 

electrode name=anode number=2 bottom 

doping uniform conc=3e 18 p.type direction=y regions=l 
doping uniform conc=3el 8 n typc direction^ regions* 3 

S 

S SECTION 3: Material Model Specification 

S 

material material=AlGaAs taupO=2.e-8 taun0=l.e-9 copt=1.5e-10 augn=5e-30augp=le-31 

impact selber anl=1.899e5 an2=1.899e5 bnl=5.75e5 bn2=5.75e5 apl=221500 ap2=221500 bp 1=657000 

bp2=657000 betan=1.82 betap=1.75 egran=0 

models material=GaAs srh auger conmob fldmob print 

$ 

$ SECTION 4: Optical source definition 

S 

beam num=l x.origin=37.5 y.origin=-1.0 angle=90.0 wavelength*. 2 min.window=-9.0 maxwindow=9.0 

S 

S SECTION 5: Initial solution 

$ 

symb newton carr=0 

solve init 

symb newton carr=2 

method trap autonr climit=75000 ctolt fact=500.0 maxtrap=10 
solve prev 

output e.field j. electron e. velocity e.mobility h.mobility qss e.temp h.temp val.band con. band qfii qfp impact 
recomb tot.doping photogen traps fiowlines ey. velocity ex.velocity hx. velocity hy. velocity u.auger u.srh u. radiative 
solve prev bl=l 

save outf=PINspec2a0.out 

S 

S SECTION 6: spectral response 

S 

log outf=PINspec2a.log 

solve prev bl=l lambda=0.2 

solve prev bl*l lambda=0.225 

solve prev bl =1 lambda =0.2 5 

solve prev bl=l lambda=0.275 

solve prev bl=l lambda=0 3 

solve prev bl=l lambda=0.325 

solve prev bl=l lambda=0.35 

solve prev bl=l lambda=0.375 

solve prev bl=l lambda =0 4 

save outf=PINspec400.out 
solve prev bl=l lambda =0.425 

solve prev bl=l lambda=0.45 

solve prev bl=l lambda=0.475 

solve prev bl=l lambda=0.5 

solve prev bl=l lambda=0.525 

solve prev bl=l lambda =0.55 

solve prev bl=l lambda =0.575 

solve prev bl=l lambda =0.6 

solve prev bl =1 lambda =0.625 


solve prev bl=l lambda =0.6325 
save outf=PIN spec6 3 2out 
solve prev bl = 1 lambda =0 65 

solve prev bl = 1 lambda =0.675 

solve prev bl=l lambda =0 7 

solve prev bl=l lambda =0.72 5 

solve prev bl = 1 lambda =0.75 

solve prev bl = 1 lambda=0.775 

solve prev bl =1 lambda=0.8 

solve prev bl=l lambda =0.825 

solve prev bl = 1 lambda=0.85 

save outf=PINspec850.out 
solve prev bl=l lambda=0.865 

solve prev bl = 1 lambda=0 875 

solve prev bl=l lambda=0,885 

solve prev bl=l lambda=0.9 

solve prev bl=l lambda=0.91 

solve prev bl=l lambda=0.9184 

save outf=PINspec2a2.out 
tonyplot PINspec2a.log 
end 



CHAPTER VI 


CONCLUSIONS 


Throughout this work, a detailed comparison of the gain and noise characteristics 
of a conventional and a doped well MQW APD was presented. The data obtained 
demonstrated a direct experimental evidence of structure induced preferential 
multiplication of electrons over holes. For the doped MQW APDs, the average gain per 
stage was calculated by comparing gain data with carrier profile measurements, and was 
found to vary from 1.03 at low bias to 1.09 near avalanche breakdown. This is in contrast 
to conventional p-i-n structures which show no gain in this regime. It was also shown that, 
as the bias was increased, the effect of the structure became less pronounced, and the 
MQW device was reduced to a conventional pin structure. Similar studies of the bias 
dependence of the excess noise characteristics show that the low-voltage gain is primarily 
due to electron ionization in the MQW APDs, and to both electron and hole ionization in 
the p-i-n APDs. Our measurements of the doped MQW APD clearly showed that for low 
gains (M < 6), the ionization ratio is greatly enhanced (k = a/p = 10 - 50) as compared to 
that in bulk GaAs (k = 1.67). At higher voltages, however, the value of k is reduced since 



the holes gain more energy from the applied electric field and are more likely to impact 
ionize. 

We’ve also observed the significant effect of mesa surface recombination processes 
on the dark current behavior of an APD. The resulting leakage currents can have dramatic 
consequences on the sensitivities and attainable gain levels in a photodiode. As a result of 
surface treatment, dark currents at low bias were reduced to as low as 1 pA. The result of 
this reduction in dark current was manifested in the structures’ high gain performance 
which exceeded 10,000 in some APDs. By being able to reduce the dark currents, it was 
possible to maintain dark current levels well below those of the photocurrents. This made 
it possible to achieve and sustain high levels of gains well beyond the onset of junction 
breakdown. 

The spectral response and quantum efficiencies for some of the structures were 
also calculated and modeled. The experimental data were very consistent with the 
theoretical models. The quantum efficiencies of the fabricated structures were relatively 
low because of the losses of photogenerated carriers due to recombination mechanisms in 
the diffusion layer. It was shown how it was theoretically possible to significantly increase 
the quantum efficiency of the devices through the introduction of a heavily doped p ++ 
GaAs top layer. Such layers help create a high-field region that will enhance the diffusion 
of photogenerated electrons toward the depletion region. 

We have also made a full investigation of the impact of doping imbalance in 
doped-well MQW APDs on the devices' IV and CV characteristics and how such an 
imbalance would affect the depletion properties of the APDs. Our theoretical models were 



in full agreement with the observed experimental data and have helped us develop a good 
understanding of some of the physical processes that take place inside a doped MQW 
APD. Our simplistic models have enabled us to easily interpret experimental IV and CV 
data and determine the extent of depletion in our APD devices. We’ve also determined 
how these parameters are affected by the p- and n-doping imbalance in the structure. The 
model predicted that a doping mismatch as small as 10% can reduce the depletion layer by 
as much as 50%. It was also shown how a large doping imbalance would cause the device 
to quickly reach avalanche breakdown in the depleted layers and would prevent full 
depletion. 

The presence of undepleted regions was also shown to be one of the major causes 
of poor time response performance in avalanche photodiodes. Partial depletion gave rise 
to diffusion-limited transient response in doped MQW structures. This was demonstrated 
experimentally to be the case by examining the change in the "diffusion tails" of the output 
pulse response of the devices as a function of applied bias. Fully depleted PIN structures 
showed fast time response even at zero bias. Doped MQW APDs showed a decrease in 
the FWHM and fall time in direct proportion to the increase in the depletion region of the 


structure. 
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APPENDIX A 


Al x Gai. x As Material System Parameters 

GaAs/AlGaAs materials are frequently used in the fabrication of heterojunction 
devices. The parameters for GaAs can, in general, be derived from those for Al x Gai- x As 
by setting the value of x to zero. In the following few sections, the various optical and 
electrical parameters relevant to GaAs/AlGaAs structures modeling will be presented. 


A.l AlGaAs Band Parameters 

The net material bandgap is chosen as the minimum value of the various default bandgap 
valleys in AlGaAs given by the following equations : 

E*(G)=Ego + x(l . 1 55 + 0.37 x) [A- 1] 

E g (L)=1.734 + x(0.574 + 0.055 x) [A- 2] 

E g (X)=l .91 1 + x(0.005 + 0.245 x) [A- 3] 

The temperature dependence of the bandgap is introduced through the Ego parameter 
calculated from the following equation 1,2 : 


E t {T) = E g (0)~ 


aT 2 

T+fi 


= £,(300) + a 


300 2 
300 + P 


T 2 

T+0 


[A- 4] 


where, 

E,(0)= 1.519 eV 
a = 5.405 x 10-4 eV/K 
P = 204.0 K 

The AlGaAs electron affinity for a given x value can be calculated from: 

X = 4.07 - 0.8 5(E g (G) - E g ( GaAs )) [A- 51 

and the density of states mass of the valence band is given by: 

m v = (0.082 3/2 + 0.45 3 ' 2 ) 2/3 [A- 6] 

A. 2 _ A l gaA s D jels ctric.Prg pgrt i es 

The default value for the static dielectric constant for AlGaAs is given as a function of x 
by: 

e’ = 13.18 + 2.9 x [A- 7] 

In the case of GaAs, x=0, and e’ = sg«a J Bo = 13.18. 

A.3 AlGaAs Bandgap Narrowing 

Heavy doping of materials usually lead to bandgap narrowing which has an important 
effect on heterojunction characteristics. According to Lundstrom 3 , the bandgap narrowing 
effects are important only for p-type regions. The default bandgap narrowing parameters 
are shown in Table 1 for AlGaAs material. 



Concentration 

cm"* 

Bandgap Narrowing 
maV 

L0E18 

31.0 

2.0E18 

36.0 

4.0E18 

44.2 

6.0E18 

48.5 

S.0E18 

51.7 

1.0E19 

54.3 

2.0E19 

61.1 

4.0E19 

64.4 

6.0E19 

61.9 

8.0E19 

56.9 

1.0E20 

53.2 

2.0E20 

18.0 


Table A-l : Default bandgap narrowing values 


A.4 AlGaAs Recombination Parameters 

The default parameters used for the various recombination mechanisms parameters are 
shown in Table A- 2 : 


Parameter 

Value 

Equation 


l.OxlO" 9 

12-191 


2.0x10-* 

(2-201 

Copt 

1.5x10’° 

[2-18J 

Cn 

5.0x10*° 

[2-17] 

Cp 

l.OxlO’ 51 

_12-171_ 


Table A- 2: Default Recombination Parameters for AlGaAs 
A. 5 GaAs Impact Ionization Coefficients 


The coefficients used in the impact ionization model described by equations [2— 21]-[2— 23] 
are provided in Table A- 3. 





























mm 

Value 

WBM 

1.899x10 s 

II 

5.75xl0 3 

A 

1.82 

mm 

2.215x10 s 

ET 

6.57x10 s 

ftp 

1.75 


Table A- 3: Impact Ionization Coefficients for GaAs 
A. 5 AlGaAs Mobility Parameters 


The default low-field electron mobility for various ranges 

of Al x Gai.*As composition is 

i 

given by the following set of equations : 



(x„ = 8000 - 1.818*10 4 x 

(0< x < 0.429) 

[A- 8] 

= 90+ 1.1435* 10 s (x-0.46) 2 

(0.429 < x < 0.46) 

[A- 9] 

p o = 90 + 3.75*10 4 (x'0 46) 2 

(0.46 <x< 0.5) 

[A- 10] 

\x« = 200-2. 0/(x-0.46) 

(0.5 <x< 1.0) 

[A- 11] 


The default concentration-dependent mobility values for GaAs are given in Table A- 4. 

The following expression is used to implement a field-dependent mobility which provides a 
smooth transition between low-field and high-field behavior. 


r 


-)«//» 


M(E) = A 


1 



1 + 


[A- 12] 














Where P is a constant with a default value of 1.0 for both electrons and holes. The 
saturation velocities for AlGaAs are given by: 

= 1.13*10 7 -1.2*10 4 7’ 
v^ = L12*10 7 -n*10 4 r 



Mobility in GaAs 
(cm 2 /v-e) 

Concentration (cm-3) 

Electrons 

Holes 

1.0*10 14 

8000.0 

390.0 

2.0* I0 1 4 

7718.0 

380.0 

4.0* I0 1 4 

7445.0 

375.0 

6.0*10 14 

7290.0 

360.0 

SO* 1 o ' 4 

7182.0 

350.0 

1.0*I0 15 

7300.0 

340.0 

2.0* I0 15 

6847.0 

335.0 

4.0»10 15 

6422.0 

320.0 

6.0*10 15 

6185.0 

315.0 

8.0*I0 15 

6023.0 

305.0 

1.0* 10“ 

5900.0 

302.0 

2.0*10“ 

5474.0 

300.0 

4.0*10“ 

5079.0 

285.0 

6.0*10“ 

4861.0 

270.0 

8.0*10“ 

4712.0 

245.0 

1.0*10 17 

4600.0 

240.0 

2.0*10 17 

3874.0 

210.0 

4.0*10 17 

3263.0 

205.0 

6.0*10 17 

2950.0 

200.0 

8.0* 10 1 7 

2747.0 

186.9 

1.0*10“ 

2600.0 

170.0 

2.0* 10“ 

2060.0 

130 0 

4.0*10“ 

1632.0 

90.0 

6.0*10“ 

1424.0 

74.5 


Table A- 4: Default concentration-dependent mobilities for GaAs 
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